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\^ ABSTRACT 

We present a joint analysis of the overlapping BLAST 250, 350, 500 |J.m, and Large 
APEX Bolometer Camera 870 [a.m observations (from the LESS survey) of the Ex- 
^•J tended Chandra Deep Field South. Out to z ~ 3, the BLAST filters sample near 

the peak wavelength of thermal far-infrared (FIR) emission from galaxies (rest-frame 
wavelengths ~ 60-200 |xm), primarily produeed by dust heated through absorption in 
star-forming clouds. However, identifying counterparts to individual BLAST peaks is 
_ ^ very challenging, given the large beams (FWHM 36 -60 arcsec). In contrast, the ground- 

k> based 870 |J.m observations have a significantly smaller 19 arcsec FWHM beam, and 

are sensitive to higher redshifts {z ~ 1-5, and potentially beyond) due to the more 
favourable negative _K'-correction. We use the LESS data, as well as deep Spitzer and 
VLA imaging, to identify 118 individual sources that produce significant emission in 
the BLAST bands. We characterize the temperatures and FIR luminosities for a sub- 
set of 69 sources which have well-measured submm SEDs and redshift measurements 
out to 2; ~ 3. For flux-limited sub-samples in each BLAST band, and a dust emis- 
sivity index (3 = 2.0, we find a median temperature T = 30K (all bands) as well as 
median redshifts: z = 1.1 (interquartile range 0.2-1.9) for ^250 > 40mJy; z = 1.3 
(interquartile range 0.6-2.1) for £'350 > 30mJy; and z = 1.6 (interquartile range 1.3- 
2.3) for 5500 > 20mJy. Taking into account the selection effects for our survey (a 
bias toward detecting lower-temperature galaxies), we find no evidence for evolution 
in the local FIR-tempcrature correlation out to z ^ 2.5. Comparing with star-forming 
galaxy SED templates, about 8% of our sample appears to exhibit significant excesses 
in the radio and/or mid-IR, consistent with those sources harbouring an AGN. Since 
our statistical approach differs from most previous studies of submm galaxies, we de- 
scribe the following techniques in two appendices: our 'matched filter' for identifying 
sources in the presence of point-source confusion; and our approach for identifying 
counterparts using likelihood ratios. This study is a direct precursor to future joint 
far-infrared/submm surveys, for which we outline a potential identification and SED 
measurement strategy. 
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1 INTRODUCTION 

Observations in the submillimetre (submm) wavelength 
band (defined here to be 200-1000 |J.m) are ideal for detect- 
ing light from massive star-forming galaxies out to cosmo- 
logical distances. It has been known since the all-sky Infrared 
Astronomical Satellite (IRAS) survey of the 1980's that such 
sources contain significant amounts of dust, so that the 
ultra-violet (UV) light of newly-formed stars is absorbed by 



the galaxies' interstellar medium (ISM) (Sanders & Mirabel 



19961. The dust is typically heated to tens of Kelvin, and 
most of the light is then thermally re-radiated at far-infrared 
(FIR) wavelengths (~60-200 (im). In the submm, the ther- 
mal spectral energy distribution (SED) drops off steeply, so 
that there is a progressively stronger negative T^-correction 
with increasing observing wavelength. The correction is so 
strong that near ~1 mm the observed flux density for a 
galaxy of fixed luminosity is approximately constant from 
l^z^m | |Blain et al.|2002| . 

Even though much of the submm band is obscured to 
ground-based observations by atmospheric water vapour, 
a number of surveys over the last decade have exploited 
transparency in several spectral windows to successfully lo- 
cate high-redshift (z > 1) dusty star- forming galaxies solely 
through their submm emission (submillimetre galaxies, or 
SMGs). Their discovery was first made with the Submillime- 



tre Common User Bolometer Array (SCUBA Holland et al 
1999[) at 850 |tm (e.g., [Small et al.|1997l [Hughes et al.|1998 



Barger et al |1998||Cowie et al |2002||Scott et al.|2002||Borys[ 
et al.[[2003[ [Webb et al.[[2003| [Coppin et al.[[2006[ ). Several 



other instruments confirmed their existence in the slightly 



more transparent 1.1-1.2 mm band (e.g., Greve et al.| 


2004 


Laurent et al.||2005| [Scott et al.[[2008| [Perera et al.[ 


2008 


Austermann et al. 


2010 


1. These ground-based surveys at 



850-1200 ixm typically cover <^ 1 deg , and detect several 
tens of sources per field. The typical angular resolution of 
these surveys is in the range ~9-20 arcsec full-width at half- 
maximum (FWHM). It is worth noting that observations at 
350 and 450 ^m have also been attempted from the ground 



(e.g.. Small et al.[[1997 Hughes et al. 


1998 Fox et al.|2002 


[Kovacs et al.||2006| [Khan et al.[[2007 


[Coppin et al.[[2008[). 



However, this wavelength range is much more difficult, due 
to increased atmospheric opacity, so that these surveys have 
only detected a handful of sources. 

While the first generation surveys successfully demon- 
strated the existence of these ultra-luminous infrared galax- 
ies (ULIRGs) at z ~ 1.5-4 (e.g|Chapman et al.|2003{[Aretx- 



aga et al. 2003 Chapman et al. 20051, sample sizes have 



been modest (typically <C 100 sources in a given field). 
Most of what is known about SMGs is based on cross- 
identifications with sources in higher-resolution data, par- 
ticularly in the radio (primarily 1.4 GHz Very Large Ar- 
Smail et al^[2000[ [Ivison et aL][2007[) and in 



ray maps, e.g.. 



the mid-IR (such as 24 (j.m Spitzer maps, e.g., Ivison et al 
[2004[ [Pope et al.|2006| l . While these counterpart identifica- 
tion strategies could be biased toward lower redshifts due 
to the positive Ti'-corrections in the radio/mid-IR, more ob- 
servationally time-consuming mm-wavelength interferomet- 



2001 Dannerbauer et al. 



ric observations (e.g., Lutz et al 
2004[[Iono et al.|2006{ [Younger et al.|2007[ ) demonstrate rea- 
sonable correspondence with proposed radio/mid-IR coun- 
terparts for a handful of sources. With accurate positions 



it is then possible to identify optical counterparts, although 
they are usually extremely faint due to obscuration by the 
same dust that makes them bright in the submm-FIR, and 
the fact that stellar light from the most distant objects gets 
red-shifted out of the optical bands into the near-IR. Obvi- 
ously, ground-based optical spectroscopy is even more chal- 
lenging given the difficulty in imaging the counterparts. 

Recent observations by the 1.8-m Balloon-borne Large 
Aperture Submillimeter Telescope (BLAST) at 250, 350, 



500 [im (a pathfinder for Herschel/ SPIKE, Pascale et al 



2008 1 toward the Extended Chandra Deep Field South 
(ECDF-S) have provided the first confusion-limited submm 
maps at these wavelengths which cover areas larger than 
Ideg^. These bands were chosen to bracket the peak 
rest-frame FIR emission from the SMG population at 
z ~ 1-4. However, given the size of its primary mirror, 
the BLAST diffraction-limited angular resolution of 36- 
60 arcsec FWHM at 250-500 \im has made associations be- 
tween submm emission peaks and individual sources at other 
wavelengths considerably more challenging than with the ex- 
isting ground-based surveys at longer wavelengths. For this 
reason many of the primary BLAST scientific results to date 
have been derived from the statistics of brightness fiuctua- 
tions for entire maps, such as the number counts (Devlin 
Patanchon et al!][2009 1 , contributions of known 



et al. 2009 



sources to the Cosmic Infrared Background (GIB Pugot ct al, 
Fixsen et aLfr998|) in the BLAST bands (Devlin ct al 



1996 



2009 



Marsden et al. 2009 Pascale et al. 20091, evolution 



in the FIR-radio correlation (Ivison et al. 20101, and the 



large-scale clustering of infrared-bright galaxies ( Viero et al, 
[2009). We emphasize that none of these results depend on 
identifying individual submm sources in the BLAST maps. 

More traditional analyses of BLAST sources identified 
through peaks in maps convolved with the point spread func- 
tion (PSF) have also been attempted ( [Dye et al.|2009[ [Dun-[ 
lop et al.|2010 Ivison et al.[[2010 l. In general, it has been a 



struggle to determine whether these peaks are produced pri- 
marily by single galaxies, or blends of several faint sources, 
necessitating either: conservative cuts in the signal-to-noise 
ratio (SNR) to consider only the very brightest sources; or 
a careful (though subjective) comparison of all the multi- 
wavelength data on a case-by-case basis to decide whether 
single or multiple objects are the likely source of the submm 
emission. In these earlier papers, Poisson chance alignment 



'P' probabilities (Downes et al. 1986| ) have been used to rank 
potential counterparts to the BLAST peaks from external 
matching catalogues, showing that 10% spurious threshold 
probabilities must be adopted to obtain reasonable source 
statistics (unlike the more conservative 5% that is typical 
in the submm community) . These methods yield limited re- 
sults for these wide-area BLAST maps, despite the fact that 
the SNR of the individual peaks rival those of most previous 
ground-based observations. 

In this paper, driven by the apparent inadequacy of 
existing methods for studying individual sources in the low- 
resolution BLAST maps, we develop improved approaches 
for: (i) filtering confused maps to find emission peaks that 
are more likely to be produced by individual (or at least a 
small number of) sources; and (ii) identifying counterparts 
to these peaks in external matching catalogues using Like- 
lihood Ratios (LR), a method which can incorporate more 
prior information than that assumed in the calculation of 
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P. In addition to deep BLAST observations of the ECDF-S, 
we also make extensive use of the deepest wide-area submm 
map at ~lmm to date: the Large APEX Bolometer Cam- 
era (LABOCA) EC DFS Submm Survey (LESS) at 870 \im 
(Weiss et al.||2009l), taken with the 12m APEX telescope 



( Glisten et al.|2006 1. This first detailed comparison between 
BLAST and longer-wavelength ground-based submm data 
helps in two key ways. First, the LABOCA beam has a 
19 arcsec FWHM (roughly half that of the 250 \im BLAST 
beam), enabling us to ascertain directly whether some of the 
BLAST peaks resolve into multiple submm counterparts. 
Second, like SCUBA, LABOCA ( |Siringo et al.|2009; ) is more 
sensitive to z > 1 sources than BLAST, and the most distant 
sources {z > 4, e.g., Coppin et al.|2009 l are expected to be 
LABOCA-detected BLAST-dropouts. Therefore this study 
will offer superior constraints on the high-redshift submm 
galaxy population than earlier BLAST studies. 

Now that we have entered the era of Herschel surveys, 
we also show that these techniques will be useful, despite 
the approximately twofold improvement in angular resolu- 
tion offered by SPIRE compared to BLAST. We explore 
this issue using simulations of SPIRE maps using the smaller 
PSFs. While we find that the situation is certainly improved 
for SPIRE, confusion will continue to seriously hamper the 
interpretation of these new surveys. 

The analysis is organized as follows. In Section [2.1| we 
summarize our treatment of the submm data using our new 
'matched filter' to identify individual peaks (full details are 
given in Appendix^. We produce an external matching cat- 
alogue in Section |2.2| combining 24 (i.m mid-IR and 1.4 GHz 
radio priors to select sources from a deep Spitzer IRAC near- 
IR catalogue. The LR identification technique is summarized 



in Section 3.1 (a full development, and calculation of priors 
are given in Appendix [B|, and it is used to produce a list of 
potential matches to the submm peaks in Section [3.2| Also, 
for cases where matches in the catalogue could not be iden- 
tified, we search for counterparts in the higher-resolution 
LESS 870 (xm peak catalogue. At this stage we have a col- 
lection of submm peaks, and a list of individual galaxies 
that we believe produce these peaks - in many cases blends 
of several galaxies. To establish their submm flux densities 
we fit PSFs at all of their locations simultaneously, in each 
of the submm maps, in Section [4.1[ The effects of confusion, 
missing identifications, and clustering are explored in Sec- 
tions |4.2||4.4| using simulations. We derive redshifts for the 
proposed counterparts in Section |5.1| and study the rest- 
frame properties of the sample in Section |5.2[ showing in 
particular how confusion may have biased some of the earlier 
BLAST results. Finally, in Section [O] we simulate SPIRE 
data to demonstrate the usefulness of our techniques for 
these new surveys. 



2 DATA 

2.1 Submillimetre Data 

2.1.1 BLAST 

In 2006 BLAST conducted a two-tiered nested survey cen- 
tred over the Great Observatories Origins Deep Surveys 
South (GOODS-S): BLAST GOODS-South Wide (BGS- 
Wide) over lOdeg^ to instrumental RMS depths 36, 31 and 



20mJy, at 250, 350 and 500 ^m, respectively; and BLAST 
GOODS-South Deep (BGS-Deep) over 0.9 deg^ to RMS 
depths 11, 9 and 6mJy. It is important to note that a sig- 
nificant additional contribution to the noise in these maps 
is produced by point source confusion, estimated to be 21, 



17 and 15mJy in the three bands (Marsden et al. 20091 



The ECDF-S is completely encompassed by the BGS-Deep 
coverage. The BLAST maps were produced using SANEPIC 
( Patanchon et al.|2008 l, and were filtered to suppress resid- 
ual noise on scales larger than approximately 10 arcmin (the 
array footprint). The BLAST beams have FWHM 36, 42, 
and 60 arcsec at 250, 350, and 500 jim. The maps and data 



reduction are discussed in detail in Devlin et al. ( 2009 1 



Details on instrument performance and calibration are pro- 



vided in Pascale et al. (20081; Truch et al. (20091 



2.1.2 LESS 

The LABOC A Survey of the Ex tended Chandra Deep Field 
South (LESS, [Wiss et al.|2009| provides deep 870 |J.m data, 
with an RMS better than 1.2 mjy across the full 30 arcmin x 
30 arcmin ECDFS. A combination of time-domain filtering 
of the raw bolometer data, as well as the suppression of 
residual noise on scales ^ 90 arcsec were incorporated as 
part of the reduction procedure. Similar to the BLAST data, 
this map was then smoothed with the 19 arcsec FWHM 
diffraction-limited PSF to identify 126 point sources in |Weiss| 
et al. ( 2009 1 above a significance of 3.7a (equivalent to a false 
detection rate of < 5%). 



2.1.3 Submm Peak Catalogues 

It is standard practice in the submm community to find 
sources in maps by cross-correlating with the PSF (this 
strategy was used both in earlier analyses of BLAST data 
and the LESS 870 vim map as noted above). This opera- 
tion is optimal for the case of an isolated point source in a 
field of statistically un-correlated noise: the cross-correlation 
gives the maximum-likelihood fiux density of a point source 
fit to every position in the map[^ However, sources are not 
isolated in the real submm maps under discussion. Their 
high surface density, combined with the large beams, prac- 
tically ensures that every pixel in these maps has at least 
some contribution from multiple overlapping sources. This 
"confusion noise", A^c, is independent of the approximately 
white instrumental noise, A^w. Confusion noise must be con- 
sidered when asking the question: what is the flux density of 
a particular point source at an arbitrary location on the sky? 
One can think of A'^c as the distribution of flux densities in a 
map with no instrumental noise source (i.e. a map of point 
sources smoothed by the PSF), precisely the distribution 



that is modelled in a P{D) analysis. As shown in Patan- 



chon et al. ( 2009 \ for BLAST, this distribution is asymmet- 



ric, with a positive tail that converges to the underlying 
differential counts distribution at large flux densities (since 
brighter sources have a lower surface density, and therefore 
stand out more against the confusion of fainter sources). For 
the remainder of this paper, both ctw and will refer to the 



^ This procedure has long been understood in astronomical data 



analysis in other wavebands, e.g.. Stetson 1 19871 
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RMS of the respective noise distributions, where the former 
is an output of the map-making process, and the latter is 
estimated from simulated maps with no instrumental noise 
added. 

The previously published BLAST PSF-smoothed maps 
lie in a regime where ctw is about a factor of two smaller 
than confusion noise, ctc. In contrast, for LESS the two noise 
components are roughly equal[^ Since confusion is a source 
of noise that is correlated on the scale of the PSF, we have 
investigated a modified filter for identifying peaks. If one 
imagines the extreme opposite case of that used to motivate 
cross-correlation with the PSF - that there is no instrumen- 
tal noise in the image whatsoever, and so the only signal is 
the confused pattern of many overlapping sources, clearly 
any additional smoothing will only make matters worse. In 
this extreme case, one must clearly bin the data into a finely- 
sampled grid and simply identify peaks without any filtering 
that might further degrade the resolution — or alternatively 
attempt to de-convolve the PSF. 

We approach this problem by developing a 'matched 
filter' that maximizes the SNR of individual point sources 
in the presence of the two noise components. The resulting 
filter in each band can be thought of as an optimal balance 
between smoothing and de-convolving the PSF, depending 
on the relative sizes of the white and confusion noise compo- 
nents. We cross-correlate the raw submm maps with these 
new filters, and produce source lists from peaks in the re- 
sulting SNR maps. A detailed description of our matched 
filter is given in Appendix with Fig. |A1| comparing the 
matched filter with the PSF at 500 vim. 

How much does this filter improve the SNR at each 
wavelength? To address this question we turn to simulations, 
drawing point source flux densities from the best-fit num- 



ber counts measured in each band (from Patanchon et al 



( 2009 1 and Weiss et al. ( 2009 1 for BLAST and LESS respec- 



tively), distributing them uniformly in maps, smoothing the 
maps with the respective PSFs, adding instrumental white 
noise to match the levels reported for the real observations, 
and finally subtracting the map means. We then compare 
the relative sizes of confusion (ctc), white (aw) and com- 
bined {at = V""! + ""w) noise components in: (i) the raw 
un-smoothed maps; (ii) maps smoothed by the PSFs; and 
(iii) maps smoothed by the matched filters. The results are 
summarized in Table [l] and in Fig. [T] a sample portion of 
the ECDF-S is compared in all four submm bands with no 
filtering, matched filtering, and PSF filtering. 

We find that the point-source sensitivities in the match- 
filtered BLAST maps improve by ~ 15-20% over those re- 



(20091 



ported for the PSF smoothed maps in Devlin et al 
(comparing the bold-face at columns for the PSF and 
Matched Filter in Table [T]). However, the increase in the 
LESS point-source sensitivities, ~ 5%, is not as significant 
due to the relatively smaller contribution of confusion. 

Finally, we produce catalogues of submm peaks for 
which we will attempt to identify the source(s) of submm 
emission. Given the marginal improvement in the 870 vim 
map using the new filter, and for the sake of simplicity, we 



For BLAST the ratios of confusion RMS to instrumental noise 
in the PSF-smoothed maps are 1.9, 1.9, and 2.4 at 250, 350, and 
500 vim. For the LESS PSF-smoothed map the ratio is 1.15. 



Table 1. Comparison of confusion (ctc), white (fw) and total 
(o"t) noise contributions in the raw submm maps, maps cross- 
correlated with the full PSF, and maps cross-correlated with the 
matched filter which compensates for confusion. Each quantity 
was estimated using simulations based on un-clustered realiza- 
tions of sources drawn from the measured number counts distri- 
butions of [Patanchon et al.| l f2009[ ) and [Weiss et al] | [2009| , for 
BLAST and LESS respectively. All noise values are standard de- 
viations in mjy. Note that our definition of confusion noise here 
is simply the RMS of a noise-free map containing point sources 
only. This calculation shows that white noise is more effectively 
suppressed by PSF filtering than matched filtering. However, the 
confusion noise is significantly larger in the PSF filtered maps 
than the match filtered maps. For this reason the total noise in 
the match filtered maps is smaller. 
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use the full 870 ^m LESS catalogue from |Weiss et al.| ( |2009 l 
(down to a SNR of 3.7a). However, we construct new 3.75a 
peak lists for BLAST (relative to instrumental noise) from 
the match-filtered maps (all submm catalogues are provided 
in Appendix [c|P] These thresholds ensure that the peaks 
are likely to be caused by submm emission rather than spu- 
rious instrumental noise; such cuts typically result in false- 
identification rates (the probability that the underlying fiux 
density within an instrument beam at that location is neg- 



ative) of order 5% (e.g., Coppin et al. 2006 Perera et al. 
|2008[ [Weiss et al. 12009] ). In total, within the region of over- 
lap, there are 64 peaks at 250 vim, 67 peaks at 350 vim, 55 
peaks at 500 vim, and 81 peaks at 870 vim. We emphasize, at 
this stage, that we do not know whether they are produced 
primarily by one, or multiple overlapping sources. We will 
will explore the properties of peaks selected this way using 
simulations in Section |4l 

2.2 Matching Catalogue 

The ideal matching catalogue for our lists of submm peaks 
would contain all of the sources that emit significantly in 
the submm, without having any spurious interlopers. The 
starting point for our catalogue is IRAC data from SIMPLE 
(gptteer IR AC/MUSYC Public Legacy in ECDF-S, [Ga^Iser] 
et al.|2006 |. This sample is approximately fiux-limited with 
Ss.(i\im <: 5/iJy and S'4.5vim <: 5 fiJy (although there is sig- 
nificant scatter at the faint end, presumably due to noise and 
completeness effects). Anecdotally, such deep near-IR cata- 
logues appear to contain counterparts to some of the faintest 
and highest-redshift submm sources detected in ground- 



based surveys (e.g.. Pope et al. 2006 Chapin et al.]|2009 



Coppin et al. 20091. However, the high surface densities, 
in this case 48.3 arcmin"^ for the ECDF-S, make it impos- 
sible to directly associate individual sources with the much 

3 The new BLAST GOODS-S match-filtered maps and peak lists 
are available at http: //blastexperlment ■ inf o/| 
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Figure 1. Comparison of maps at all four wavelengths (rows), and smoothing scales in a O.Sdeg X 0.3 deg patch of the ECDF-S. The 
first column shows un-smoothed maps (noisy, but diffraction-limited resolution), the second column maps smoothed by the matched 
filter (greatly improved SNR at the expense of a slight degradation in the resolution), and the last column maps smoothed by the PSF 
(less improvement in the SNR, and resolution degraded by \/2). The greyscales indicate the significance of point-source flux densities in 
the maps, ranging from — Scr (white) to +4o- (black), considering both the instrumental and confusion noise contributions to each pixel 
(estimated from simulations — see boldface values in Table [ij . 
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lower-resolution submm peaks. For this reason it is generally 
necessary to use other priors to cull known spurious sources 
(such as stars, or galaxies with little or no dust). Fortunately, 
in this field there also exist deep 24 \im Spitzer and 1.4 GHz 
VLA maps. These two wavelength regimes have been used 
considerably in the past to identify submm sources as they 
are sensitive to the presence of dust and star-formation ac- 
tivity, respectively. We check for emission in these additional 
data sets to reduce the 33962 sources from SIMPLE cover- 
ing the 703 arcmin^ region of overlap (Fig. [2]| , to a total of 
9216 entries, a surface density of 13.1 arcmin"^. There are 
8833 sources detected at 24 ^m, 1659 at 1.4 GHz, and 1276 
at both 24(i.m and 1.4 GHz. 

In Sections [2m and l2T2l we describe our treatment of 
the mid-IR and radio data, respectively. Then, in addition 
to the previous plausibility arguments, we use stacking to 
check directly whether our combined catalogue reproduces 
the diffuse measurements of the GIB, and hence determine 
how complete it is in the four submm bands under consid- 
eration in Section [2.2.31 



2.2.1 FIDEL 

A deep Multiband Imaging Photometer for Spitzer (MIPS) 
24 and 70 [im catalogue was produced by [Magnelli et al.| 
(20091), combining Spitzer data from GOODS-S with the 



Far-Infrared Deep Extragalactic Legacy survey (FIDEL, P.I. 
Mark Dickinson). The PSFs have FWHM 5.5arcsec and 
16 arcsec, in each band respectively. The catalogue uses SIM- 
PLE as a positional prior enabling de-blending of sources 
down to separations as small as 0.5 times the MIPS FWHM. 
This catalogue is the same that was used for other recent 



BLAST studies in EGDF-S pevlin et al.||2009| [Dye et al. 
2009| [Marsden et"aL|p009| [Pascale et al.||2009[ ), and is con- 



sistent with producing the entire GIB across the BLAST 
bands. This catalogue is also used in the LESS identifica- 
tion paper (Biggs et al. submitted), although it is known 
that such catalogues do not reproduce the entire GIB at 
wavelengths approaching ~ 1mm (Wang et al. 2006 Pope 



2007 Marsden et al. 2009 1 



2.2.2 Raa 



The VLA 1 .4 GHz data are from the survey of jMillerj 
|et al.| ( |2008[ ). Given the low elevation of this southern 
field, and observing from the VLA in the north, the syn- 
thesized beam is significantly elongated with dimensions 
2.8 arcsec x 1.7 arcsec. We use a new radio catalogue pro- 
duced by Biggs et al. (submitted) that was developed to 
identify counterparts to LESS sources. An initial catalogue 
of 3 (T peaks in the map is produced, and then Gaussians are 
fit at each of those positions, allowing the sizes to vary. Since 
this fit is particularly noisy at the faint end, we cull sources 
with integrated flux densities that are less significant than 
2 a. Finally, we only include sources that lie within 2 arcsec 
of the IRAG positions. This strategy enables us to go signif- 
icantly deeper than the published 7 a catalogue fr om jMillerj 
et al. (20081, or the 5 a catalogue from Dye et al. (20091, at 



the expense of missing a handful of brighter radio sources 
that do not appear to have IRAG associations. Since Monte 
Carlo simulations based on the submm data and matching 



catalogue are used in Section [3. 2j to calculate the probabil- 
ity that individual counterparts are real, any spurious ra- 
dio sources near the detection threshold will simply reduce 
the identification efficiency (see also Fig. B3l. Finally, we 



note that Biggs et al. (submitted) also find identifications 
for LESS sources using a combination of the radio, MIPS, 
and IRAG catalogues, although the IRAG data are not ex- 
plicitly used as a prior for the radio positions. 



2.2.3 Stacks 



We use the method of Marsden et al. (20091 to measure 



the submm surface brightnesses of the complete SIMPLE 
catalogue, as well as the 24 ^m (FIDEL), 1.4 GHz (VLA), 
and combined 24 ^.m and 1.4 GHz (FIDEL-I-VLA) subsets 
of SIMPLE in Table [2] These values are compared to the 
absolute measurements of jFixsen et al.|fl998| ). 

First, we find that the stack on the complete (and very 
high-surface density) SIMPLE catalogue results in lower val- 
ues than the lower-surface density FIDEL and FIDEL-KVLA 
catalogues. Even though we know that this IRAG catalogue 
contains sources that are not strong submm emitters (e.g., 
following the arguments in Section 3.1 of Marsden 



stars) 



et al. ( 2009 1 we would not expect these extraneous sources 



to bias the result unless they were somehow correlated on the 
sky (they are simply an additional source of noise). There- 
fore, we believe that this peculiar result simply demonstrates 
the presence of an un-characterized systematic in the SIM- 
PLE spatial distribution (i.e., they are anti-correlated with 
submm emission). We note that Marsden et al. (20091 also 



found unexpected behaviour when stacking on a high-surface 
density optical catalogue. 

Second, we find that stacks on the FIDEL catalogue 
yield results consistent with those quoted in jMarsden et alT] 
( 2009 1 (our values are slightly lower, as expected, because we 
have not corrected for completeness in FIDEL at faint flux 
densities as they did). Within the uncertainties, this cata- 
logue is consistent with reproducing the entire GIB at 250- 
500 [im, but recovers less than half of the GIB at 870 vim, 
consistent with the previous stacking measurements noted 
earlier. 

Finally, we see that the addition of 1.4 GHz detected, 
but 24 (j,m un-detected sources ("VLA only"), to FIDEL sys- 
tematically increases the value of the stack slightly; i.e. the 
stack on FIDEL-I-VLA is greater than the stack on FIDEL 
by itself. We use this fact as our primary justification for 
including the additional faint radio sources. 

While this catalogue arguably contains the majority 
of the sources that produce significant submm flux in the 
BLAST bands, the catalogue is certainly missing a signifl- 
cant portion of the 870 \im emitters. Unfortunately, we do 
not know whether these missing sources are, on average, 
the same sources that produce the 870 [im peaks, or fainter 
sources that do not typically contribute to the brighter 
peaks. We note that the LESS identification paper. Biggs et 
al. submitted, also fails to identify counterparts for a num- 
ber of the 870 |J.m sources using similar matching data. We 
will attempt to address the impact of this shortcoming in 
Section [43l 
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Figure 2. Relative coverage of data sets in ECDF-S. The background greyscale image shows the BLAST 250 ixm match-filtered SNR 
map (this field is completely encompassed by the BGS-deep region described in [Devlin et al.||2009[ l, scaled between — lOcr (white) and 
+13(7 (black). The solid white contours show the LESS SNR map at levels 3, 6 and 10 it. The white dashed lines show the LESS 1.3 and 
2.2 mjy instrumental noise contours. The dashed black line indicates the 708arcmin'^ region common to the VLA and FIDEL survey 
coverage within which we perform out counterpart search. The ECDF-S presently has the best (widest and deepest) submm coverage 
from 250-870 ixm on the sky, with the mid/far-IR and radio data of matching quality required to identify counterparts. 



3 CROSS-IDENTIFICATIONS 

Given the poor positional uncertainties inherent to the cur- 
rent generation of submiUimetre waveband surveys (typi- 
cally several arcsec), there are usually many potential op- 
tical counterparts for each SMG. Thus it has usually been 
necessary to search for identifications in lower surface den- 
sity catalogues at radio and mid-IR wavelengths. Also, as 
mentioned in Section [2. 2| emission in these two wavebands 
are expected to be physically correlated with the submm- 
FIR emission; there is not a similarly strong correlation 
for optically-selected sources. The method usually adopted 
is to estimate 'P' Poisson chance alignment probabilities 



( Downes et al.|1986 l in order to exclude the least likely can- 
didates (although this does not provide a probability that 
a given source 'is the counterpart'). This calculation only 
uses the source counts of the matching catalogue and an 
empirically derived maximum search radius. The expected 
distribution of offsets for true matches is not used, except 
perhaps to set the search radius. 

In this paper we take a different approach, using a 'like- 
lihood ratio' (LR) formalism. The basic idea attempts to an- 
swer the following question: given a potential counterpart to 
the submm peak, what is the relative likelihood that it could 
be a real counterpart given its measured properties (e.g., ra- 
dial offset, flux density, colour etc.), versus the probability 
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Table 2. Surface brightnesses re sulting from stacks on different catalogues as compared to absolute measurements of the total CIB 
(values shown in brackets are fromj Fixsen et al.[l998| , in order of decreasing surface density (quoted values in arcmin"^): SIMPLE is the 
entire IRAC catalogue; FIDEL+VLA includes sources from SIMPLE that exhibit 24 |j.m emission, and additional VLA 1.4 GHz sources 
that have a significance > 2a within 2 arcsec of a SIMPLE position; FIDEL is the 24 i^m catalogue based on SIMPLE from |MagnelIi] 
|et al.| | |2009l l; and VLA is the subset of SIMPLE sources with 1.4 GHz emission that are not members of FIDEL. The units of the surface 
brightness measurements are nWm~^sr~^. The FIDEL+VLA catalogue (indicated in boldface) is used for matching throughout this 
paper. 



Catalogue Surface 250 |J.m 350 |xm 500 [im 870 |^m 

Density (10.4 ±2.3) (5.4 ±1.6) (2.4 ±0.6) (0.47 ±0.1) 



SIMPLE 


48.3 


4.33 


± 


0.87 


3.71 


± 


0.51 


1.75±0.30 


0.13 


± 


0.01 


FIDEL+VLA 


13.1 


8.25 


± 


0.46 


4.60 


± 


0.27 


2.09 ±0.16 


0.17 


± 


0.01 


FIDEL 


12.6 


8.16 


± 


0.45 


4.56 


± 


0.26 


2.02 ±0.16 


0.17 


± 


0.01 


VLA only 


0.54 


0.35 


± 


0.10 


0.16 


± 


0.06 


0.15 ±0.03 


0.004 


± 


0.001 



that it is a chance interloper (given the background source 
counts as a function of measured properties)? Versions of 
this technique have been used in a variety of contexts (e.g., 
Sutherland fc Saunders]p92l |Mann et all[T997l |Rutledge| 



et aL 2000 1 . Clearly the LR can explicitly use information 
such as the expected positional uncertainties, whereas the P 
calculation does not (although the maximum search radius 
implicitly incorporates some of this information). It should 
be noted that, in the past, colour-based priors have been 
used to cull matching catalogues, and hence reduce the sur- 
face density of spurious sources, before finding counterparts 



using P statistics (e.g Pope et al. 2006 Yun et al. 2008 



Biggs et al. submitted). By itself the LR can still only be 
used to rank potential counterparts (similar to the problem 
with P statistics) , but we attempt to establish both the false 
identification rates, and identification completeness rates for 
given absolute LR thresholds in each band. 

It appears that the reason the LR formalism has not 
been used for submm identification work in the past is due 
to its reliance on prior information which historically has 
been extremely difficult to estimate (see discussion in 'Ser-' 

With individual 



jeant et al. 2003 



Clements et al. 



20041 



surveys covering ^Ideg and having typically fewer than 
several tens of peaks per field, and very low SNR, the pre- 
cise positional uncertainty distribution is unknown. How- 
ever, the BLAST and LESS data, combined with the deep 
radio and Spitzer mid-lR data covering the ECDF-S, enable 
estimates of priors (such as the radial offset distribution of 
counterparts) with sufficient precision to produce useful re- 
sults. Forthcoming Herschel, SCUBA-2, and LMT surveys 
will have better angular resolution, depth, and cover sub- 
stantially larger areas, so that the methods employed here 
will also be fruitful (however, note that these future surveys 
will depend on radio and mid/near-IR data of comparable 
area and depth for identifying counterparts, and they do not 
presently exist). In Section 3.1 we summarize the LR method 
and the priors that we have developed. In Section [3^ we then 
use this method to identify a list of potential counterparts 
from our matching catalogue to the submm peaks. 



3.1 Likelihood Ratios and priors 

As described in Section |2] we have produced a matching 
catalogue based on sources from SIMPLE (IRAC) that ex- 
hibit either mid-infrared [Spitzer MIPS photometry from 



FIDEL), or VLA 1.4 GHz emission. For the purpose of iden- 
tifying counterparts we have focussed on three features of 
this catalogue in addition to positions: 24 [j.m and 1.4 GHz 
fiux densities (when available) — both of which have com- 
monly been used in the calculation of P values; and the 
[3.6] — [4.5] IRAC colour, which is sensitive to redshift 



(e.g.,|Simpson & Eisenhardt|1999 Sawicki 


2002 Pope et al. 


|2006| |Yun et al.||2008| |Devlin et al.||2009 


as it traces the 



peak of the rest-frame stellar bump. 

We fully develop the LR formulae and priors in Ap- 
pendix |B] but provide the main results here. Given the flux 
densities 5* (at 24 (i.m and 1.4 GHz in this case), the IRAC 
colour c, and the distance r to the jt\i matching catalogue 
source from the ith submm peak for which we are searching 
for a counterpart, we calculate the LR: 



Li, 3 — 



(1) 



Here a characterizes a radially-symmetric Gaussian posi- 
tional uncertainty, q{S, c) is the prior distribution for flux 
densities and colours of matches to submm peaks, and p(S, c) 
is the background source distribution. 

All of the priors, cr, q{S,c) and p{S,c), have been es- 
timated directly from the data by counting sources in the 
matching catalogue as a function of each property around 
submm peak positions (of order ~60 in each submm band), 
and comparing to the counts for the entire matching cat- 
alogue over the full survey area. We find that sources in 
the matching catalogue near submm peaks (i.e., potential 
counterparts): have radial offset distributions that are pro- 
portional to the instrumental PSF sizes, as expected; tend to 
have 24 yuai and 1.4 GHz flux densities that are brighter than 
typical sources in the catalogue; and have redder [3.6] — [4.6] 
IRAC colours than typical sources in the catalogue, par- 
ticularly at 500 and 870 \im. We also find that, on average, 
there are multiple extra sources from the matching catalogue 
near each submm peak, an excess E, of: 3.2 at 250 (xm; 3.4 
at 350 ^m; 3.7 at 500 ^.m; and 2.2 at 870 ^tm. This result 
demonstrates that the submm data are highly confused, but 
clustering in the matching catalogue may also have an im- 
pact. 

There could be submm-faint matching catalogue mem- 
bers that cluster around a smaller number of submm-bright 
sources which produce the observed submm peaks. However, 
we note that E is strongly correlated with the PSF size; if 
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Table 3. Summary of the matching catalogue selection function 
and priors for Likelihood Ratios. Flux cuts for the matching cata- 
logue are only approximate (see Section |2.2| . No hard cutoffs are 
used as priors, there are only weights applied as a function of the 
distance between the submm peak and proposed counterparts, 
the 24 i^m and 1.4 GHz flux densities, and [3.6] — [4.5] near-IR 
colours (see Appendix [b] and Figs. [B1]B4| 



Catalogue selection 


Notes 


524|xm <; 13/iJy 


Spitzer/MWS, FIDEL 




(uses SIMPLE IRAC positions as a 




prior, see below) 




or 


•S'i.4GHz ^ 20/iJy 


VLA 




and 


Ss-evim ^ 5/iJy 


Spitzer/IRAC, SIMPLE 


'S'4.5|J.m ^ S^tJy 





Prior 


fir) 


Favour nearby counterparts 


'j('S'24n.m) 


Favour brighter S24[im 


9('5l.4GHz) 


Favour brighter S'i.4ghz 


9(c) 


Favour redder [3.6] — [4.5] colour 



this clustering scenario were the dominant effect we would 
expect the same integrated excess regardless of beam size. 
We will further explore the potential impact of clustering in 
Section Ol 

We believe that our excess counting procedure yields a 
high-SNR measurement of the positional offset distribution, 
since we need only bin measurements for the submm peaks 
(in each band) along one coordinate, r, which we then fit 
with a simple one-parameter model, /(r) = {r/a^)e~^ ^^'^'^ . 
However, there is an implicit assumption that the sources 
(both in the submm and matching catalogues) are spatially 
un-clustered. In Table [3] we summarize the matching cata- 
logue selection function, (described in Section 2.2 I, as well 
as the basic effect of the priors used in the LR calculation. 

Ideally we would bin q{S, c) and p{S, c) along all 
three axes simultaneously (two flux densities and a colour), 
but in practice this is not feasible given the numbers of 
submm peaks that we have to work with. We therefore 
handle the priors independently, e.g., estimating q{S,c) ~ 
q{S2i)q{Sr)g{c), and p(5, c) ~ p{S24)p{Sr)p(c). This as- 
sumption certainly introduces a bias, since in practice these 



properties are correlated (for example, see Fig. 7 in Dye 
|et al.| |2009[ showing the correlation between 24 fim and 
1.4 GHz flux densities for counterparts to bright BLAST 
sources). These correlations are simply a reflection of the 
fact that submm galaxies have a particular range of SED 
shapes (including the radio and near-IR discussed here). We 
have not attempted to measure these shapes, but an alterna- 
tive method might consider a range of plausible, physically- 
motivated SEDs as part of the identiflcation process (e.g. 



Roseboom et al. 2009) . Instead, we have chosen to com- 



pensate for the bias by using Monte Carlo simulations to 
estimate a threshold in the LR that produces a false identi- 
flcation rate of 10%. 

We establish the appropriate level by choosing 10,000 



random positions in the fleld, and calculating the LR of 
matching catalogue sources around those positions. In each 
submm band we choose a threshold in LR that rejects 90% 
of the matched sources from this random sample. For conve- 
nience, we then calculate normalized LRs for counterparts 
to submm peaks in each band by dividing the raw LR from 
Eq.[l]by these thresholds. We therefore only consider poten- 
tial counterparts to submm peaks those objects for which the 
normalized LR is greater than 1. While this choice of nor- 
malization is arbitrary, the relative LRs for potential coun- 
terparts are meaningful (i.e., a LR of 2 indicates double the 
relative likelihood that it is real compared to a LR of 1). We 
also estimate that a proposed counterpart with a normal- 
ized LR of 1 is about three times more likely to be real than 
spurious. 

Finally, note that the LR only gives the relative chance 
that a proposed counterpart is real. We have not attempted 
to derive an absolute reliability, 7?, the probability that a 
proposed counterpart is correct, as in [Sutherland fc Saun-| 



ders (19921, since their formulation assumes that there is a 
single counterpart to each peak (see Appendix [B| . Instead 
we rely on our threshold in LR to ensure our spurious frac- 
tion of 10% for the ensemble of proposed counterparts. 

In summary, while our simplifying assumption that the 
flux densities and colour are independent of one another is 
incorrect, we have established a cut on LR that will restrict 
the number of false positive identiflcations in the matching 
catalogue to around 10% (with respect to the number of 
submm peaks). This idea of setting a threshold to reject un- 
related sources is similar to adopting a cut on P, although 
here we have included more prior information to improve the 
efficiency. However, we warn the reader that the calibration 
of the LR and P are both tied to our assumption that clus- 
tering in the matching catalogue around submm positions is 
a negligible effect. 

3.2 Potential counterparts 

We search for counterparts to all of the submm peaks inde- 
pendently in all four submm bands out to a search radius of 
60 arcsec. In practice we could use a much larger search ra- 
dius, but, since the expected radial probability density, /(r), 
rolls off to ~ before this radius, even at 500 \im, there is no 
difference in the list of counterparts with normalized LR > 1 
if it is increased. We identify the following numbers of poten- 
tial counterparts in each band: 52 at 250 \im; 50 at 350 vim; 
31 at 500 [im; and 66 at 850 vim. 

As noted in the previous section, we expect to en- 
counter several counterparts, on average, for each submm 
peak. However, with the threshold LR chosen, we only de- 
tect a fraction of the total number of counterparts expected: 
22% at 250 ^m, 19% at 350 ^m, 13% at 500 |i.m, and 32% at 
870 [im. It is probably the case that this low identification 
rate is simply due to the quality of the submm data. 

For comparison, we have also identified potential coun- 
terparts using P statistics based on the 24 [im and 1.4 GHz 
fiux densities. We counted the numbers of IDs with P < 0.1 
using both of these sub-catalogues around random positions 
(as a control, similar to the method employed to normalize 
the LR), and around submm positions (similar to Fig. 3 in 
[Chapin et al.|2009| ). We discovered that, while the radio cat- 
alogue was well-behaved (a 10% false ID rate was obtained 
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for random positions), we obtained too many false IDs with 
the 24 |i.m catalogue. This result is probably demonstrat- 
ing that the 24 yun catalogue is slightly clustered. We there- 
fore tuned the cut on P24 to 0.08 to obtain the desired 10% 
rate using random positions. The search radii that we used 
were 1.5 (the single parameter in the radial offset distri- 
butions); these search radii are roughly comparable to those 



adopted in previous studies of BLAST peaks (Dye et al 



2009 Ivison et al. 2010 Dunlop et al. 20101, and encom 



pass approximately 68% of the true counterparts inferred 
from the excess counting statistics. Using P24 < 0.08 we 
would find 45, 31, 12 and 33 potential counterparts, and 
using Pr < 0.10 we would find 51, 35, 17 and 44, at 250, 
350, 500, and 870 vim, respectively. Clearly there is a signifi- 
cant (although modest) improvement using our LR calcula- 
tion over P statistics, particularly at the longer wavelengths 
where the [3.6] — [4.5] IRAC colour is a good discriminator 
of redshift. 

As our goal is to study the properties of BLAST- 
selected peaks, we identify all of the unique sources from our 
matching catalogue that potentially produce the observed 
250, 350 and 500 (xm emission. However, we only consider 
matches to the LESS peaks that lie within a. 2 a search ra- 
dius of any 3.75a BLAST peak (combining both the BLAST 
and LABOCA positional uncertainties in quadrature). With 
this search radius we expect to find 95% of the real matches, 
and given the surface density of the LESS catalogue, we will 
have a 10% spurious ID rate (the same as that adopted for 
the cut on LR). In total, 42 out of the 81 LESS peaks that 
land within the survey region are associated with BLAST 
sources, and 36 of them are identified in the matching cat- 
alogue using LRs. The remaining 6 LESS peaks are still 
included as potential matches for BLAST peaks, although 
they may themselves also be blends of multiple galaxies 
within the LABOCA beam, and none of them have asso- 
ciated radio, mid- or near-IR flux densities with which to 
conduct further analysis. 

The end result of our matching procedure is a list of 
118 unique sources that are believed to contribute to the 
submm peaks in all four bands (including the 6 that are 
simply LESS 870 vim sources with no matches at other wave- 
lengths). The coordinates of these sources, and the submm 
peaks to which they were matched, are given in Table |C5| 
Postage stamps showing the locations of the matches in re- 
lation to the submm positions are shown in Fig. |D1[ For 
each ID, in addition to the LR, we also provide P values for 
matches to the 24 [xm and 1.4 GHz radio catalogues, using 



search radii 1.5ar (see Table Bll 



Note that the third columns in the lists of submm peaks. 
Tables |Clf|C4[ give references to the individual identifica- 



tions from the matching catalogue in Table C5 (again, within 



a 2(T search radius). Since these sources may have been 
matched to submm peaks in any band, the total numbers 
of counterparts listed here exceed the numbers of matches 
found independently in each band. In other words, we can 
now see how the simultaneous observations in the four 
submm bands have helped one another: 77%, 72%, 76%, and 
86% of the 250, 350, 500, and 870 [im peaks have at least 
one potential counterpart identified. In many cases the same 
sources appear in several bands, and the higher-resolution 
observations have enabled counterpart identifications where 
the lower-resolution observations failed. However, counting 



the total number of proposed counterparts, we find 64, 65, 
59, and 51 sources at 250, 350, 500, and 870 These 
numbers still fall far short of the total numbers expected, 
especially in the BLAST bands. We will attempt to quan- 
tify the impact of the missing matches in Section [4. 3[ 



4 SUBMILLIMETRE SPECTRAL ENERGY 
DISTRIBUTIONS 

In this section we re-measure the submm flux densities from 
the combined list (i.e., selected in any BLAST band) at the 
positions of their proposed counterparts in the raw submm 
maps. Using simulations, we explore the impact of confusion, 
missing sources, and clustering on these measurements. Fi- 
nally, we fit these observed-frame SEDs with simple isother- 
mal models, and measure the observed number counts in our 
catalogue as a function of limiting fiux density. Together, 
these calculations allow us to explore bias and completeness 
effects for our sample. 



4.1 Submm photometry 

Under the assumption that the potential counterparts iden- 
tified in the previous section produce the observed submm 
emission (and are not simply other galaxies clustered around 
the submm peaks), we return to the four submm maps and 
perform simultaneous fits of point sources at all 118 loca- 
tions to measure their fiux densities. This procedure is ex- 
pected to reduce the Eddington-like bias, or flux-boosting, 



inherent to low-SNR submm surveys ( Coppin et al. 2005 \ 



in two ways. First, since peaks are initially selected in three 
different bands, the component of bias introduced by instru- 
mental noise is reduced (as it is independent in each map) 
— submm peaks preferentially detected on positive noise 
excursions in one map will not necessarily also land on pos- 
itive noise excursions in other maps. Second, since we allow 
for the possibility of multiple counterparts to each submm 
peak (a hypothesis confirmed by the excess counting statis- 
tics described in Section Bll, the simultaneous fit can, to 
some extent, de- blend some of the brighter confused sources 
(confusion itself also contributes to Eddington bias, and un- 
like the instrumental noise, is correlated between the submm 
bands). 

The fit is performed by modeling the emission of the 
counterparts as the submm PSFs scaled by their unknown 
submm flux densities Si at the locations from the matching 
catalogue. Under the assumption that the instrumental noise 
in our submm maps is un-correlated from one map pixel to 
the next (a reasonable assumption for the raw, un-smoothed 
maps on the scale of the PSF) , there is a simple maximum- 
likelihood solution for the Si that takes into account the cor- 
relations that arise in cases where multiple sources overlap 
within a PSF footprint — we follow the derivation provided 
in Appendix A of Scott et al. (20021. This solution only 



^ Note that at 870 |J.m 57 matches are indicated, but 6 of those are 
the LESS sources themselves, leaving 51. We have only searched 
for counterparts to the 42 LESS peaks that appear to be associ- 
ated with BLAST peaks, so the average counterparts per peak is 
actually 51/42=1.21 
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uses the maps, instrumental noise estimates, and source po- 
sitions. No preferential weight is given to counterparts with 
larger LRs, so the noise in the answer only depends on how 
well the map is fit using our simple parameterization. 

A downfall of this approach is that we ignore the ad- 
ditional component of confusion noise from un-identified 
sources that is correlated on the scale of the PSF. We there- 
fore estimate the total noise by adding the confusion noise 
for the simulated raw maps from Table [T] to the variances 
for each source flux density, al — Cov{Si, Si). This opera- 
tion should give good estimates for isolated sources, but we 
warn that it produces an under-estimate of the variances for 
the most confused sources. In Section 14.31 we will test the 
validity of our estimated uncertainties using simulated data 
sets. 

For isolated sources, the recovered flux density is identi- 
cal to that obtained from the PSF-smoothed map at the lo- 
cation of its counterpart, and its value is un-correlated with 
the measured flux densities for all other sources in the map. 
However, for blended sources, the total flux density in the 
map is divided among the multiple counterparts, and there 
are non-negligible covariances Cov(S'i, Sj) for all sources i,j 
that lie roughly within a FWHM of each other. We therefore 
evaluate the full expression for the covariances between mea- 
sured flux densities, i.e., using the off-diagonals of Eq. All 
inlr •■'^^ 



Scott et al. (20021. 



The individual observed-frame submm SEDs based on 
these measurements are given in Table \C6\ We also display 
the submm SEDs in Fig. |El] along with the 1.4 GHz, MIPS 
24 and 70\im flux densities (when available), and the pho- 
tometry from the 3.6, 4.5, 5.8 and 8.0 iJ.m IRAC bands (all 
sources) from the matching catalogue. 

4.2 Confusion 

Our method assumes that the proposed counterparts to the 
submm peaks comprise all of the galaxies that contribute 
significant submm emission. However, we have made two 
fairly arbitrary choices: we select only peaks that have a sig- 
nificance of 3.75 a over the instrumental noise levels; and we 
only consider sources with a 10% threshold false association 
rate from the LR analysis. To assess the impact these choices 
have on the measured flux densities and completeness, we 
have generated simulated maps drawing sources from the 
measured number counts in the BLAST bands from lPatan-l 
|chon et al.| ( |2009[ ), and then adding appropriately scaled 
white noise to mimic the estimated instrumental noise lev- 
els. 

First, we identify individual peaks in the simulated 
BLAST maps above the same 3.75cr SNR threshold as for 
the real data using the same matched filter. Given the sizes 
of the BLAST beams, and the surface density of the sources, 
every location in the filtered maps has a contribution from 
multiple submm galaxies (even if they are extremely faint). 
For each 3.75cr peak we therefore identify the single source 
that makes the largest contribution to the observed flux den- 
sity from the input catalogue at that location — consider- 
ing the PSFs in each band, input source brightnesses and 
their distances from the peaks in the filtered maps. In this 
way a faint source will only be identified provided that it is 
very close to the submm peak in question, and exceeds the 
brightnesses of the tails of all the more distant sources in 



the catalogue. We then re-run this procedure on 100 inde- 
pendent realizations of the maps at each wavelength to fully 
characterize the scatter in the results. Each time we ran- 
domly select a different differential counts distribution from 
the actual Markov Chains produced from the the P{D) fits 
p009|. 
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We find that these brightest sources statistically con- 
tribute fractions 0.73t!5 j^, 0.57lg:^^, and 0.59,t°o%l (means 
and 95% confidence intervals) of the peak flux densities in 
the filtered maps at 250, 350 and 500 \im. Note that the frac- 
tion can be greater than one since the simulation is noisy 
(the source may have landed on a large positive noise ex- 
cursion). This test shows us that, using a 3.75cr cut, submm 
peaks are usually a significant blend of two or more individ- 
ual sources, although there is an incredibly large scatter; a 
peak may have many contributors, but it is also true, in some 
cases, that a peak is dominated by a single bright source. 
This result is broadly consistent with a similar set of simula- 



tions used by Moncelsi et al. (20101 to correct BLAST peak 



flux density biases. Also, we have noted that the radial dis- 
tribution of the brightest sources identified for each submm 
peak (not shown) broadly resemble the radial distributions 
estimated for the LR analysis (Eq. B2|. 



There are numerous obvious examples of blends of 
sources from the matching catalogue in the real submm 
maps. We have flagged 42/118 of the most extreme cases 
with the letter 'C in Fig. |E1| indicating that they are con- 
fused to the point that the submm photometry cannot be 
used reliably, particularly in the BLAST channels. For ex- 
ample, sources 2, 3, 4 and 5 comprise one of the most con- 
fused regions of BLAST emission in the entire ECDF-S, as 
can clearly be seen in the postage stamps (Fig. Dll. By 



comparison, the superior LESS resolution can nearly resolve 
the entire feature into a string of individual peaks. The in- 
ferred flux densities at 250-500 \uxi therefore have strong 
anti- correlations, since the emission from those four sources 
must sum to the total integrated flux density of the feature. 
In fact, the maximum-likelihood solution we have adopted 
can even allow negative values, a problem which occurs in a 
number of the most confused examples. The low SNR of the 
submm maps, combined with confusion from fainter submm 
sources, and the close proximity of the counterparts, has 
resulted in flux densities with drastically under-estimated 
error bars in cases such as sources 2-5. In contrast, sources 
61 and 62 are an example where the joint-fit at the po- 
sitions of two nearby counterparts has recovered plausible 
flux densities in all the submm bands (this is a low-redshift 
interacting pair first discussed in Dunlop et al.|2010 l. In this 
case, the BLAST SNR is much higher, and the two potential 
counterparts have sufficient separation to disentangle them. 

Next, we investigate the completeness by counting the 
number of sources above different flux limits in the input 
catalogue that are recovered in the match-filtered source list 
(again, considering only the single brightest submm galaxies 
that contribute to the observed brightness). The recovered 
percentages are: 50% above 30 mjy, and 90% above 60 mjy 
at 250 [J.m; 50% above 15mJy, and 90% above 45mJy at 
350 [im; and 50% above 10 mJy, and 90% above 25 mJy at 
500 nm. 
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4.3 Missing counterparts 

Another significant problem that we face is the issue of miss- 
ing source matches to the submm peaks. As mentioned in 
Section |2.2.3| the matching catalogue is probably missing 
a significant fraction of the 870 \im emitters, and possibly 
a smaller fraction of the 250-500 \im emitters. More impor- 
tantly, in Section |3.2| we were not able to identify a large 
portion of the counterparts to the submm peaks expected in 
the matching catalogue. Therefore, the measured flux densi- 
ties for the sources that were correctly identified are noisier, 
and perhaps biased, due to these missing sources in the fit- 
ting procedure. We attempted to account for this noise in 
Section |4.1| by adding the RMS of simulated, instrumental 
noise-free maps in quadrature to the noise returned from the 
fitting procedure. Here we will use simulations to estimate 
how biased and noisy this procedure is. 

We use the simulated maps from the previous section 
that were generated using realizations of sources drawn from 
the measured counts distributions for the real maps, and 
assigning them random positions. We now also simulate 
870 nm maps using the best-fit counts as reported in |Weiss| 
et al. ( 2009 1. For each realization we add a random 20% un- 



certainty to the total number of sources drawn from the dis- 
tribution to approximate the error indicated in the faintest 
bin of their cumulative catalogue-based counts. 

To approximate the source identification procedure, we 
first produce 3.75(t peak lists, and then identify the brightest 
sources from the input catalogue that contributed to each of 
the peaks. We include as matches only those sources from 
the input catalogue that contribute more than a threshold 
fraction of the observed peak, chosen to produce the same 
average number of matches per peak as we obtained using 
the real maps and matching catalogue. In this way we as- 
sociate a range of matches from the input catalogue (with 
known flux densities) to each observed peak, with a similar 
surface density of matches as for the real data. Ifowever, we 
stress that this is only a plausible simulation, since there is 
no guarantee that the matches proposed for the real data are 
in fact the brightest contributors as is the case for this simu- 
lation. A significantly more complicated simulation could be 
undertaken in which we: generate galaxies with full radio- 
submm-IR SEDs (that are consistent with the true surface 
densities of sources in each band); create a matching cat- 
alogue (i.e., simulating radio, mid- and near-IR catalogues 
with realistic noise); repeat the process of estimating priors; 
and finally use LRs to propose matches for each peak. How- 
ever, we felt that such a simulation was beyond the scope of 
this paper, and opted instead for the simpler approach that 
captures most of the necessary ingredients. 

With simulated maps, and lists of proposed matches 
to each peak, we repeat the maximum-likelihood fitting op- 
eration of Section |4.1[ in all of the submm bands. Since we 
know the true flux densities for each of the matched sources, 
St , we are able to directly probe the scatter of the observed 
flux densities, R = (St — So)/o'c where So is the inferred 
flux density, and CTo its uncertainty derived from the fitting 
process (which accounts for instrumental noise and overlap 
with other nearby sources), and then adding the additional 
confusion noise for the simulated raw maps from Table [T] m 
quadrature. The expectation is that this distribution has a 
mean of zero, and a standard deviation of one. Combining 



the results for all 100 simulations in each band, we measure 
the following mean values and standard deviations for R: 
0.17 ± 1.23 at 250 Hm; 0.19 ± 1.20 at 350 nm; 0.25 ± 1.20 at 
500 ^m; and -0.05 ± 0.97 at 870 \im. 

These results suggest that there is only a small up- 
ward bias in the measured flux densities across the BLAST 
bands, of order a"o/4, and negligible bias at 870 \im. In abso- 
lute terms, we can scale this bias to approximate flux den- 
sity units by multiplying by the mean values of (To in each 
band: -1-3.2 mjy at 250 jxm; -1-3.2 mjy at 350 vim; -f3.5mjy 
at 500 vim; and — 0.05mJy at 870 [im. We also flnd that our 
inclusion of the confusion noise only under-estimates the 
noise in the BLAST bands by at most ~23% (as discussed 
in Section [4.11 our estimate is a lower-limit on the noise for 
any particular source). 

Finally, as noted at the end of Section [2]23j our match- 
ing catalogue is potentially very incomplete at 870 \im. How- 
ever, the results of our 870 [im simulation suggest that the 
missing sources are so faint that they do not contribute sig- 
nificantly to the submm peaks we are analyzing. Again, for 
this to be true, the sources that are in the matching cata- 
logue need to account for the bulk of the brightest 870 ^im 
emitters in the sky. 

4.4 The impact of clustering 

In addition to the confusion arising from a uniformly dis- 
tributed population of submm emitters (i.e., chance super- 



positions of objects at different redshifts, as in Section 4.2 1, 
the submm emitters themselves, and/or the matching cata- 
logue, could also be clustered. There are three distinct cases 
that would affect the results in this paper worth considering: 

(i) The low-resolution submm peak could be resolved 
into multiple components at approximately the same red- 
shift. This is plausible, since it is known that about 10% 
of SCUBA sources are associated with double radio sources 



(e.g., 


Ivison et al. 


|2002| Chapman et al. 2005 


2006 


Ivison et al. 


|2007|), and Vaisiinen et al. 



Pope et al. 



20101 show 



that many of their 180 (i.m-selected sources are blends of 
multiple galaxies at 2; < 0.3. We also note that a signifi- 
cant clustering signal on angular scales < 1 arcmin has been 
measured for the LESS catalogue in excess of the Poisson 
expectation ( Weiss et al.|2009 l. 

(ii) The submm peaks could instead be dominated pri- 
marily by single matching catalogue sources with a lower 
surface density, in which case the extra sources counted in 
the radial excess plots may be spatially correlated, but are 
not otherwise directly associated with the submm emitter 
(e.g., galaxies with lower star-formation rates in the same 
structures) . 

(iii) Massive foreground structures could enhance the 
brightnesses of background galaxies through lensing. This 
scenario would have a similar effect to the previous one; 
there could be a number of additional foreground galaxies 
near the positions of submm peaks, even though they do not 
themselves contribute significantly to the submm flux. 

The impact of these clustering scenarios is not easy to 
assess accurately with simulations because it is simply un- 
known how all galaxy populations cluster throughout the 
history of the Universe; indeed, this is one of the major 
outstanding questions in modern Cosmology. While there 
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are a number of ways to simulate clustered galaxy popula- 
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tions, such as the Halo model (Mo & White||1996l used to 



fit BLAST data in Viero et al. f2009), or more complicated 



semi-analytical models (e.g. |Baugh et al. 2005J , it is beyond 
the scope of this paper to test the full range of models that 
are plausible. On the other hand, there is clear evidence that 
some clustering is required to explain both the BLAST and 
LESS submm maps. Since we have argued that the bulk of 
the submm emission is produced by sources in the match- 
ing catalogue, a reasonable approach is to incorporate these 
real, and clustered positions in our simulated maps to see 
what impact they have. 

We repeat the simulations of the previous sections, 
drawing sources from the measured counts distributions, but 
now randomly assigning them positions from the real match- 
ing catalogue. With this limitation, we can only draw the 
same number of sources as in the real matching catalogue, 
so we restrict ourselves to a flux-limited sample that results 
in the same surface density. We then produce maps, and 
identify 3.75ct peaks as before. 

To assess the relative impact of clustering, we also pro- 
duce simulations with the same surface density of input 
sources, but using uniformly distributed positions. In both 
cases we measure the radial excess counts in the match- 
ing catalogue around peaks, with respect to the entire cata- 



logue (as in Figure Bl|. We find that the clustered simula- 
tions yield systematically larger excesses per submm peak, 
by factors 1.07, 1.02, 1.02, and 1.24 at 250, 350, 500, and 
870 \im, respectively. The trend is for this excess to be larger 
for the highest-resolution measurements, which might be 
expected given the large 2-point correlation measurements 



from Weiss et al. ( 2009 1 growing towards scales < 1 arcmin. 



However, considering the scatter in the 100 simulations, we 
also find that the uncertainties in the excess measurements 
are 21%, 14%, 19%, and 16% in the four bands. In other 
words, when testing the hypothesis that the difference ob- 
served for the clustered simulations are significant compared 
to the un-clustered simulations, only the excess at 870 vim 
is marginally significant (a 1.5cr outlier). 

However, we warn that our simulation is only testing 
a particularly weak form of clustering. It is probable that 
subsets of our matching catalogue have significantly dif- 
ferent angular clustering signatures when compared to the 
catalogue as a whole, and it may be possible to identify 
them through their redshifts, colours, and brightnesses. If 
such populations are correlated more, or less strongly with 
the submm emitters, their may be additional significant bi- 
ases in our LR approach - effectively CEises (ii) and (iii) de- 
scribed above. We note that Chary & Pope (20101 investi- 



gate the impact of just such an effect on stacking analyses in 
this field, although the results depend substantially on their 
model for the redshift evolution of submm galaxies and their 
multi-wavelength SEDs. 

We note additional support for the hypothesis that clus- 
tering has a negligible impact in Section |3.2| the excess 
matching catalogue counts around submm positions are cor- 
related with the angular resolution of the submm maps. The 
same integrated excess would be measured, regardless of the 
PSF shape, if the submm emission were produced by single 
galaxies. 

As a final word on this subject, what would the effect 
be on our analysis if there were significant clustering in the 
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Figure 3. The integral number counts for our sample (solid lines), 
compared with the total counts inferred from the P(D) analysis of 
[Patanchon et al.| ( |2009| (dotted lines), at 250 (blue), 350 (green) 
and 500 |xm (red) . 



matching catalogue around the submm peaks? In this case 



the radial excess distributions (Fig. Bll, and our measure- 
ments of q[S, c) (Figs. B2 ■ B4 1 would trace the properties of 



the clustered objects (e.g., spatial extent, colours and bright- 
nesses), rather than the properties of submm sources, and we 
would end up with many more false positives than the target 
10% rate. While we do not believe this is the most likely sce- 
nario based on the arguments made in this section, and the 
fact that the SEDs for the counterparts seem to match mod- 
els for star-forming galaxies, only higher-resolution studies 
will be able to settle this issue unambiguously. 



4.5 Isothermal SED models and number counts 

We fit optically-thin isothermal modified blackbody func- 
tions (modified blackbodies henceforth), 5*1^ oc v^B^{T), to 
the new submm photometry, and use Monte Carlo simula- 
tions to characterize the uncertainties (as described in Sec- 
tion 4 of Chapin et al.||2008 l. Following the work of Wiebe 



et al. (20091, who examined the detailed spatially-resolved 



SEDs of several nearby resolved galaxies observed with 



BLAST, as well as an earlier study by Klaas et al. (20011 



who combined ground-based submm photometry with FIR 
measurements of ULIRGs, we choose to model the submm 
emission using an emissivity index of /3 = 2.0. The only 
free parameters are the amplitudes and temperatures (see 
Table [C6] and Fig. |E1| . 



To test the completeness of our catalogue, and also to 
gauge the degree to which our procedure has dealt with fiux 
boosting in the submm peak catalogues, we compare the 
integral source counts from our sample (summing the num- 
ber of sources in the catalogue above a given flux density 
limit, and dividing by the survey area of 708 arcmin^) to 



the total counts inferred from P(D) analysis in Patanchon 



et al. (20091. Rather than using the submm photometry di- 



rectly, we use the fiux densities from the fitted SED models 
evaluated in each band. The results of this comparison are 
shown in Fig. [S] In the 250 and 350 [im channels the cata- 
logue counts slightly exceed the P{D) counts above approx- 
imately 30 and 40 mJy, respectively. While this excess shows 
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there is still some influence from boosting, its effect has been 
drastically reduced when compared with the individual flux- 



limited BLAST catalogues (see Fig. 11 in Patanchon et al. 



2009). Below these levels the catalogue is clearly incomplete 
at 250 and 350 ^.m as the counts rapidly flatten. At 500 ^.m 
the catalogue counts have a similar qualitative shape, but lie 
below the P{D) counts at all flux densities; the complete- 
ness is about 66% above 20mJy. This result is expected 
given the more limited success we have had in identifying 
counterparts at 500 \im (see Table Bl I. Note that these ap- 
proximate completeness estimates are consistent with the 
simulations described at the end of Section 14.21 



5 DISCUSSION 
5.1 Redshifts 

Many of the proposed counterparts have either optical spec- 
troscopic or photometric redshifts in previously published 



catalogues 


(Wolf et al. 


2004 


2008 Grazian et al. 


2006 


Brammer et al. 2008, Rowan- Robinson et al.||2008| |Taylor| 


et al.||2009 


\. For those counterparts that do not, the IRAC 



colours may be used as a crude redshift estimator. We use 



the redshift catalogue from Pascale et al. ( 2009 1 which com- 



bines the various photometric and spectroscopic redshfts 
in the literature with the BLAST redshift survey of |Eales| 
et al. (20091, and then we add additional redshifts identi- 
fied in more recent BLAST follow-up studies ( [Ivison et al.| 
2010[|Dunlop et al.|2010l [Casey et al.|2010[ ). In Fig.|El[ the 
redshifts are indicated with 's', 'p' or 'i', indicating spec- 
troscopic, optical photometric, or IRAC-based photomet- 
ric redshift measurements, respectively. In total, there are 
76/118 sources with usable submm photometry. Of those, 
69 have redshift estimates: 23 are optical spectroscopic red- 
shifts; 35 are optical photometric redshifts; and 11 are IRAC 
photometric redshifts. The full list of redshifts that we have 
adopted is given in Table [CH] We warn that the IRAC-based 
photometric redshifts are highly uncertain on an object- 
by-object basis, and are biased low for the higher-redshift 
{z > 2) sources (see Fig. 4 in [Pascale et al.|2009p . 

In Fig. |4|we show the redshift distribution for these 69 
sources, as well as sub-samples using flux density-limits in 
each of the BLAST bands corresponding roughly to the flux 
densities at which the counts begin to turn over significantly 
in Fig. [3] (as a rough proxy for the point at which complete- 
ness begins to drop). The median of the entire distribution 
is z = 1.1 with an interquartile range 0.3-1.9. Note that if 
we exclude the the IRAC-based photometric redshifts the 
median of the entire sample increases slightly to z = 1.3. 

This redshift distribution is qualitatively similar to the 



deep 250 |am survey of Dunlop et al. (20101 in GOODS-S at 



the center of the ECDF-S, but shows a significantly greater 
tail of sources beyond 2 = 2 compared to the shallower sur- 



vey of Dye et al. ( 2009 1. This latter discrepancy is probably 
due to a combination of increased depth in our submm cat- 
alogue, and better completeness in the high-redshift coun- 
terpart identifications. 

The fiux-limited distributions clearly show a trend from 
low to high redshift with increasing wavelength in Fig. [4] 
a median z = 1.1 with an interquartile range 0.2-1.9 at 
250 [i.m; a median z = 1.3 with an interquartile range 0.6- 
2.1 at 350 |J.m; and a median z — 1.6 with an interquartile 
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Figure 4. The redshift distribution for the 69 non-confused 
sources with redshift estimates. Also shown are the redshift dis- 
tributions for fiux-limited sub-samples, 5250 > 40mjy, S350 > 
30 mjy, and 5500 > 20 mjy, chosen to correspond approximately 
to where the counts in our catalogues flatten significantly com- 
pared to the total population - sec Fig. [3| The asterisks indicate 
the medians of these sub-samples. 



range 1.3-2.3 at 500 (xm. This trend is consistent with the 



results from BLAST stacking analyses (Devlin et al. 2009 



Marsden et al. 2009 Pascale et al. 2009 1 which show that the 



CIB is produced by higher-redshift galaxies with increasing 
wavelength. 

With redshift estimates for our proposed counterparts, 
we are also able to check whether clusters of sources con- 
tributing to single submm peaks lie at a single redshift, or 
at a range of redshifts. If, as we have asserted, most of our 
submm peaks are chance superpositions of un-related galax- 
ies, we would expect most clumps of proposed counterparts 
to lie at different redshifts. On the other hand, clusters of 
proposed counterparts at the same redshift would be consis- 
tent with the clustering scenarios described in Section [4. 4| 

There are cases of what appear to be groups of pro- 
posed counterparts to single peaks at different redshifts. See, 
for example, sources 2-5 in Fig. |D1[ In this particular case 
there also appears to be a fifth significant source of emis- 
sion in most of the submm maps that is not identified, this 
being west and slightly south of the main clump, perhaps 
coincident with a faint radio source that lies within the sat- 
urated source apparent in the 3.6 \im map. The clump of 
sources 13-15 is a similar example. Sources 52 and 53 form 
a more well-separated example, with the former peaking in 
the 250 vim map, the latter in the 500 [im map, and with 
a double peak in the 350 ii.m map. The two optical photo- 
metric redshifts appear to be significantly different (0.1 and 
0.6, respectively). Some of these sources could also be ex- 
amples of foreground galaxies lensing background sources, 
although this effect is more difficult to quantify without ac- 
curate spectroscopic redshifts for most of the sample, and 
lensing models for each case. 

There are also examples of peaks that could plausibly 
be interacting pairs at the same redshift (case (i) from Sec- 
tion 4.2 1, more along the lines of radio-doubles detected in 



SCUBA surveys. The clearest example is the low-redshift 



interacting pair of sources 61 and 62 (also noted in Dunlop 



BLAST and LABOCA observations of the ECDF-S 15 



et al. 112010 1. Some more possible examples are: 74 and 75; 



76 and 77; 83 and 85; and 113 and 114. 



5.2 Rest-frame SEDs 

The SEDs for all 118 sources are shown in Fig. |E1[ For the 
69 sources with redshift estimates and useful photometry, we 
convert the observed temperatures from the modified black- 
body fits of Section|4]to rest-frame temperatures, allowing us 
to probe the cold dust SEDs of the sample. The distributions 
of these temperatures for the entire sample, and flux-limited 
sub-samples in each BLAST band are shown in Fig. [5] The 
total distribution has a median T = 29K and interquartile 
range 23-36 K. For reference, re-running the analysis with 
P — 1.5 increases the temperatures by about 5K. 

There is almost no variation in the temperatures of 
the different flux-limited sub-samples: a median T = 30K 
with an interquartile range 25-39 K at 250 |i.m; a median 
T = 30 K with an interquartile range 23-35 K at 350 vim; 
and a median T = 30K with an interquartile range 22- 
38 K at 500 vim. The reason for this can be seen in the 
bottom panel of Fig. [5] even though the longer-wavelength 
channels tend to be biased to selecting higher-redshift (and 
hence more luminous and warmer) galaxies, the extra galax- 
ies picked up at low redshift (and hence lower luminosity) 
by the short-wavelength channels are biased toward warmer 
temperatures in those volumes. We have also checked for bi- 
ases in the total temperature distribution due to the IRAC- 
based photometric redshifts, but do not find a significant 
difference for the z > 1.5 galaxy temperatures when sources 
with these photometric redshifts are excluded. 

We infer cold-dust temperatures that are generally 
warmer than the mean temperature of 23.4 K and 1 a width 
of 5.2 K (also assuming /? = 2.0) reported for the robust, but 
shallower. 



BLAST sample in |Dye et al.| ( |2009| ). While this 



discrepancy may partly be due to the relative depths of the 
samples, there is also strong evidence that our treatment of 
confusion and the inclusion of LABOCA 870 vim data have 
resulted in less-biased temperatures. 

The 870 vim data place a firm constraint on the 
Rayleigh- Jeans part of the SED. Since the longer-wavelength 
BLAST data are more confused (and hence more bi- 
ased), there is a systematic preference for shallower spec- 
tra, or cooler apparent temperatures. However, the higher- 
resolution LESS data are less prone to this problem, and 
therefore greatly improve the fits (if a smooth SED model 
such as the modified blackbody is assumed). Sources 6, 23, 
25, 34, 36, 44, 48, 67, 104, and 118 in Fig. [ET] are more 
extreme examples that show this effect clearly. In all these 
cases, the SED model is a reasonable fit to the 250, 350, 
and 870 vim data points, but the 500 vim data lie above the 
model. 

Restricting ourselves to the the 26 galaxies at z ^ 1.5 
with usable submm photometry we find a median temper- 
ature of 36 K. This is in good agreement with Fig. 9 from 
Dunlop et al. (20101 which shows that their 525o/5870 data 
are consistent with a /3 = 1.5, T = 40K modified blackbody 
at these same redshifts (remembering that if we had chosen 
P = 1.5 instead of 2.0 our temperatures would be larger by 
about 5K). 

We also fit the star-forming galaxy SED templates of 
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Figure 5. The dust temperature distribution, assuming Sv oc 
u^''^ Bi,(T), for the 69 non-confused sources with redshift esti- 
mates. The bottom panel shows redshifts plotted against temper- 
ature, and the top panel collapses the redshift axis into a his- 
togram. Similar to Fig. |4] we indicate the distributions for flux- 
limited sub-samples in each BLAST band, and indicate the me- 
dians with asterisks. We indicate the sources with IRAC-based 
photometric redshifts, as they are highly uncertain, and may be 
biased low. If these sources are at higher redshifts, they will move 
up and right in the lower panel. There is a negligible trend in 
rest-frame temperature for sources selected in different BLAST 
bands. 



70 Vim) . Since these provide a fairly restricted range of SED 
shapes, they are only used to illustrate several basic features 
in the radio and mid/near-IR SEDs of our sample with re- 
spect to the rest-frame FIR peak (i.e., they should not be 
considered accurate flts beyond the 870-70 vim data). First, 
they give an indication of the location of the ~ 1.6 vim stel- 
lar bump that is redshifted into the IRAC bands for many of 
our sources (i.e., if the wavelength of the bump in the model 
corresponds to a local maximum in the IRAC data the red- 
shift estimate is plausible). Also, since these are templates 
for star-forming galaxies, sources from our sample with large 
excesses in either the mid-IR or radio may harbour AGN. 
Of the 69 sources with usable submm photometry and red- 
shifts, 6 clear examples (8%) seem to exhibit such excesses: 
source 6 (mid-IR excess); source 34 (radio excess); 49 (radio 
and mid-IR excess); 51 (mid-IR excess); and 52 (radio ex- 
cess); and 67 (mid-IR excess). For the remaining 42 submm 
confused sources, we fit the Dale et al. ( 2001 ) SED templates 



Dale et al. (2001 1 to the submm and FIR photometry (870- 



to all of non-confused data (VLA and Spitzer) to illustrate 
what plausible values in the submm might be. Finally, we 
note that some of the sources with IRAC-based photomet- 
ric redshifts have particularly poor SED fits in the near-IR 
when compared to these templates, such as sources 33, 46, 
50, and 117, as well as particularly cool inferred rest-frame 
temperatures (< 20 K). It is quite likely that these photo- 
metric redshifts are biased low, resulting in under-estimates 
of the inferred temperatures, and luminosities. 
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Next, we normalize each of the full radio-submm-IR 
SEDs by their total IR (TIR) 10-1000 ^.m luminosities us- 
ing the modified blackbody fits. The sample is then plotted 
in bins of Ltir in Fig. |6] This strategy minimizes scatter 
in the region of the SEDs probed by the submm photome- 
try, and enables us to examine trends in the relative radio 
and mid/near-IR scatter as a function of luminosity. For 
reference, the same TIR-normalized SED from |Dale et al.| 
(20011 (with logjQ(5'6o/'S'ioo) = —0.2) is also plotted over 
each luminosity range to highlight differences. In the first 
instance this figure demonstrates that the vast majority of 
our proposed multi-wavelength counterparts are plausible: 
while there are a handful of significant outliers, the relative 
intensities of the radio and IR emission to the rest-frame FIR 
peak show reasonable consistency with each other, and with 
the models. It is clear that the radio-FIR correlation holds 
for most of the sample since the reference FIR-normalized 
star-forming galaxy SED has a radio spectrum that passes 
through most of the VLA data in each luminosity bin, re- 
gardless of the wavelength of the FIR peak (consistent with 
the findings of Ivison et al. '2010'). On the other hand, the 
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mid/near-IR SEDs exhibit significantly larger scatter; par- 
ticularly in the range 11.5 < log(L) < 12.8 (a factor of ~ 100 
peak-to-peak). Most of the sources mentioned previously 
with mid-IR excesses land in the two luminosity bins that 
span this range. While a correlation between AGN activity 



and FIR luminosity has been observed before (e.g. Takeuchi 



et al. 2003 2004 Valiante et al. 20091, it is curious that 



our most luminous 12.8 < log(L) < 13.5 bin exhibits a more 
compact spread in SED properties consistent with pure star- 
formation. This result may not be significant, however, since 
this last luminosity bin has only 10 sources (whereas the 
next two fainter bins have 26 and 20 sources respectively). 
Fig. [5] also shows a strong correlation between luminosity 
and temperature: the least luminous sources peak at wave- 
lengths A > 100 |J.m, while the most luminous sources peak 
at wavelengths A < 100 \im. Finally, this plot shows that the 
more luminous sources have a greater ratio of FIR luminos- 
ity to starlight (the 1.6 (xm bump); the lowest-luminosity 
sources have SEDs more closely resembling normal star- 
forming galaxies than ULIRGs. 

To explore the correlation between luminosity and tem- 
perature more fully, and the potential for its redshift evolu- 
tion, we compare our sample with the local Universe. We 
use as our reference the distribution in integrated 42.5- 
122.5 [im FIR luminosity, and C = logj^Q(S'6o/>S'ioo) colour 



as measured by Chapin et al. (20091 based on IRAS data 
(an updated version of the analysis in Chapman et al.|2003 1. 
At redshift ~ 2.5 the entire BLAST 250-500 ^m bandpass 
closely matches the 42.5-122.5 (im coverage from the IRAS 
60 and 100 [im channels (see Fig. 6 in Chapin et al.|[2009 l, 
so that tests for evolution at that redshift have only a min- 
imal dependence on the SED model used to fit the data. 
However, at lower redshifts, the BLAST bandpass samples 
significantly longer wavelengths in the rest-frame. Since the 
modified blackbody models used to fit the submm data 
fall-off much faster in the mid/far-IR than for real galax- 
ies (since the shorter wavelengths sample warmer and/or 
optically-thick dust), we have fit the submm and FIR data 
for our sample using the more realistic star-forming galaxy 
SED templates of [Dale et ah] ( |2001[ ). In Fig. [7|we plot the 
colours and luminosities for our sample using these SED 
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Figure 6. SEDs of the 69 galaxies with non-confused submm 
photometry and redshift estimates, normalized by the total in- 
frared (TIR) 10-1000 |xm luminosities of their modified blackbody 
fits. The SEDs are plotted in bins of Ltir and offset by multi- 
ples of 103 for clarity. The SEDs of individual galaxies are shown 
in Fig. |E1| A single LxiR-normalized SED template from |Dale| 
|et al.| ( |2001^ , for which logj^Q(S6o/5ioo) = —0.2, is shown as a 
dotted black line for reference in each bin. This comparison high- 
lights several features: more luminous sources peak at shorter 
wavelengths (they are warmer); the radio-FIR correlation (the 
normalized model passes through the radio data irrespective of 
the FIR peak wavelength); the spread in the near/mid-IR SEDs 
(largest in the bins spanning 11.5 < log(L) < 12.8); and an in- 
crease in the size of the FIR peak with respect to stellar light (the 
1.6 |xm bump) with increasing luminosity. 



fits. We colour-code the sources by redshift bin: z < 0.5 
(blue crosses); 0.5 < z < 1.5 (grey asterisks); 1.5 < z < 2.0 
(green triangles); and 2 < z < 42 (red squares). Note that 
there is some quantization along the vertical axis for the 
sample since there are only 64 SED templates in the li- 
brary, with colours ranging from about —0.55 to 0.26. Since 
this range is not quite as broad as that observed for the 
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real sources, there is also some clipping (mostly problematic 
for the coolest sources at the bottom of the plot which are 
simply assigned the template SED with the most negative 
colour). We then compare this distribution with the local- 
Universe measurement of Chapin et al. ( 2009 1 , including the 



1(7 envelope (solid and dotted black lines respectively). 

We find that the temperatures for our sample lie sys- 
tematically toward cooler (more negative) values than in the 
local Universe. At low redshifts, the lower temperatures are 
expected, since the submm bands sample the Rayleigh- Jeans 
tail of the thermal emission and are therefore biased toward 
the detection of cooler sources for a given flux-density-limit. 
Similarly, observations on the shorter- wavelength side of the 
thermal emission peak (such as with IRAS) are biased to- 
ward the detection of warmer sources — an effect which was 
accounted for in [Chapin et al.| ( |2009[ ) in order to infer the 
properties of the entire population. Selection effects alone 
may therefore explain much of the discrepancy in temper- 
atures seen here, and for 850 jim (e.g., [Kovacs et aL]|2006[ 
Coppin et al.|2008 1 selected galaxy populations (see also the 
discussion in Section 4.3 of Dye et al.|2009 1. 

To test this theory, we have used a simple simulation 
to explore the effects of our selection function on the tem- 
peratures and luminosities of galaxies in different redshift 
slices. We take the local FIR colour-luminosity distribution, 



clearly reproduced. This result suggests that there is no sig- 
nificant evolution in the correlation between FIR luminosity 
and temperature at high redshifts, at least for the most lu- 
minous (I/FiR <; 10^^ Lq) galaxies probed by our sample at 
z ^ 1. 

While the lack of evolution in our sample appears to 
be at odds with the evolution observed in the stacked SEDs 
with redshift described in Pascale et al. ( 2009 1, we note that 



our sample contains only ULIRGs at 2 > 1.5, whereas the 
stacked SEDs are dominated by significantly fainter galaxies. 
It is possible that 'normal' star-forming galaxies are gener- 
ally warmer in the early Universe than in the present-day, 
while massive starburst galaxies do not differ appreciably. 
Another possible explanation is simply that the complete- 
ness to less luminous (and presumably cooler) galaxies drops 
off faster at high redshift than it does for the more luminous 
(and warmer) galaxies in Pascale et al. ( 2009 1 — resulting 
in a bias to warmer stacked temperatures. 

Our finding may also have interesting consequences for 
recent studies which indicate that SMGs are generally more 
extended and cooler than local ULIRGs. These claims have 
been made based on a variety of observations, including re- 



solved MERLIN radio morphologies (e.g. Chapman et al. 
|2004|, Spitzer MIPS and IRAC SED s (e.g. [Hainline et al 



2009| , and Spitzer IRS spectra (e.g. Menendez-Delmestre 



$(L,C), from|Chapin et al.|([2009| and assume that it does |et al.||2009[ ). Such observations point to a systematic dif- 



not evolve with redshift. We then calculate the total num- 
ber of objects as a function of L and C in a redshift slice 
zo < z < zi as J^^ ^{L,C){dV/dz)dz, where dV/dz is the 
differential volume element, and then use the [Dale et al.| 
( |2001[ ) SED templates to extrapolate to the BLAST wave- 
lengths, and hence estimate observed flux densities. While 
the relative numbers of galaxies in different redshift slices 
will be incorrect, since it is known that significant luminos- 
ity and/or density evolution in the local luminosity func- 
tion is required to fit observed submm number counts (e.g.. 



ference in the physical conditions of major-mergers which 
produce SMGs (e.g., see resolved molecular gas observations 
and discussion in Tacconi et al.|2008 \ , compared to ULIRGs 



in the local Universe. We note that Chanial et al. (20071 



Rowan-Robinson|2001| Scott et al.|2002| Lagache et aL|2003| 



Valiante et al. 2009| Rowan-Robinson |2009| Wilman et al.| 



2010), this calculation should give a good idea of the selec- 



tion function within a redshift slice (provided that evolu- 
tion across the slice is negligible). We then apply flux den- 
sity cuts of 40, 30 and 20mJy at 250, 350, and 500 [im re- 
spectively to mimic our observational selection function. We 
show 90% confidence intervals for these model distributions 
in the z < 0.5 and 2.0 < z < 4.0 redshift bins with solid blue 
and red contours in Fig. [7] showing that they extend signif- 
icantly below the rest-frame colour-luminosity correlation. 
To illustrate how this selection happens, we show the effect 
of the BLAST flux limits in the colour-luminosity plane at 
a single fixed redshift z = 0.5 as a thick dashed line (sources 
at that redshift should only be observable to the right of the 
limit). The detected galaxies tend to pile-up against this 
limit, since there are many more galaxies at fainter lumi- 
nosities than at brighter luminosities in a given volume. If 
this limit were a vertical line in the plot, we would then ex- 
pect the observed luminosities and colours to cluster around 
the rest-frame correlation (solid and dotted lines) . However, 
since the limit is inclined, a dis-proportionate number of 
objects are detected toward the lower- left (lower luminosi- 
ties, and cooler temperatures). Despite the simplicity of the 
model, the large noise in the observed sample, and poten- 
tial for incorrect redshifts, the general trends in the data are 



found a positive FIR luminosity-size correlation with a sim- 
ilar form to the temperature-luminosity correlation that we 
have used here. It may therefore be the case that the low- 
temperature bias for submm samples has picked out galaxies 
that are also more extended. 



5.3 Predictions for Herschel/SPIRE surveys 

Will the techniques described in this paper be useful for 
the new generation of Herschel/ SPIKE surveys? To ad- 
dress this question, we have again turned to simulations, 
creating realizations of sources with flux densities drawn 
from the BLAST P{D) counts. These sources are assigned 
uniformly-distributed random positions covering the same 
area as the ECDF-S observations described in this paper, 
and then smoothed with the smaller 18, 25, and 36arcsec 
FWHM Gaussians to approximate the SPIRE beams. We 
have chosen a pixel size of 6arcsec, and added Gaussian 
noise with RMS 1.8, 2.3, and 5.0 mjy at 250, 350, and 
500 [im, respectively. We find that the RMS due purely to 
point source confusion is 7.0, 6.8, and 6.6 mJy at 250, 350, 
and 500 \un respectively. This simulation is a reasonable 
approximation of the SPIRE Science Demonstration Phase 
observations of GOODS-N, taken as part of the Herschel 
Multi-tiered Extra-galactic Survey (HerMES). Initial num- 
ber counts were shown to be consistent with the BLAST 
results in I Oliver et al. ( 2010 1, and our instrumental and con- 



fusion noise estimates are close to those reported in [Nguyen] 
et al. (20101. A 0.1 deg x 0.1 deg subset of the simulated 



250 [im data are shown in the left panel of Fig. |8] 

Next, we produce match-filtered maps, following the 
prescription of Appendix |A] The total noise (combining 
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Figure 7. Rest-frame logj^Q(S'6o/Sioo) colours vs. integrated 42.5— 122.5 |J.m far-infrared (FIR) luminosities inferred from fits of ' Dale| 
|et al.| | |2001[ l SED templates to the submm and FIR photometry for our sample. The solid black line shows the local correlation derived 
from IRAS galaxies in [Chapin et al!] ( |2009^ , with the dotted black line its 68% confidence interval. The symbols are data points from 
our survey, divided into sources at 2: < 0.5, 0.5 < z < 1.5, 1.5 < 2 < 2.0, and 2.0 < z < 4.02 (including the highest-redshift source 
in the sample). We indicate the sources with IRAC-based photometric redshifts as they are highly uncertain, and may be biased low. 
If these sources are at higher redshifts, they will move up and right in this plot. For comparison, we generated a realization of sources 
from a non-evolving FIR colour-luminosity distribution l |Chapin et al.| |2009), extrapolating to the BLAST wavelengths using the [Dale] 
|et al.| | [2001[ l SED templates, and applying flux density cuts 40, 30, and 20mjy at 250, 350, and 500 [xm respectively. Sub-samples were 
then made for the 0.0 < z < 0.5 and 2.0 < z < 4.0 redshifts bins, and their 90% confidence intervals are shown as solid blue and red 
contours respectively. For further illustration we also show the BLAST selection function at z = 0.5 with a thick black dashed line. The 
reasonable agreement between the model contours and the sample demonstrates that the generally cooler (more negative) colours for 
our sample can be explained by selection effects, rather than evolution in the colour-luminosity correlation at high rcdshift. 



white and confusion noise) is estimated to be 5.8, 5.6, and 
5.7 mjy at 250, 350, and 500 mjy, respectively. We com- 
pare the 250 Jim match-filtered map with the raw and PSF- 
filtered maps in Fig. [S] In all three cases we indicate submm 
peaks with brightnesses 5*250 > 10 mJy (which are > 5 cr 
detections with respect to instrumental noise), as well as 
the locations of the known input sources to the simulation. 
Clearly these data have a SNR (with respect to instrumen- 
tal noise) significantly higher than the BLAST 250 vim data; 
therefore the improvement in the match-filtered results over 
the raw and PSF-filtered maps is more pronounced than for 
the BLAST data shown in Fig. [l] 

We draw attention to the fact that the SPIRE peaks 
are obviously blends of multiple significant sources in most 
cases. To quantify the confusion we have repeated the radial 
excess counts following Appendix |B1| using the known in- 
put source positions as the matching catalogue. At 250 vim, 
which we emphasize is the least confused SPIRE band, we 
find: ar — 4.0arcsec (radial positional uncertainty) and 



E = 1.8 (average number of excess sources contributing to 
peak) in the raw maps; ar = 2.3arcsec and E = 1.3 in 
the match-filtered maps; and Or ~ 5.4arcsec and E — 2.3 
in the PSF-filtered maps. In other words, peaks in the raw 
250 [im maps are often blends of about two galaxies, while 
the matched-filter is able to significantly de-convolve the 
beam and improve the positional uncertainties. This exam- 
ple also clearly demonstrates that convolution with the full 
PSF is the wrong approach to take with deep SPIRE data. 

We have repeated this procedure for the longer wave- 
lengths (also using a limiting fiux density of lOmJy for 
the peak catalogues), and find for the match-filtered maps: 
ar = 4.1arcsec and E — 1.8 at 350 (im; and ar = 6.8arcsec 
and E = 2.5 for the match-filtered maps at 500 \mi. There- 
fore, confusion is a significant issue, even in the match- 
filtered maps, at the longer wavelengths that are better 
suited to detecting the highest-redshift sources. 

In summary, deep SPIRE data are considerably con- 
fused at levels that lie between the LABOCA and BLAST 
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Figure 8. A simulated 0.1 dcgx 0.1 deg 250 \mi SPIRE observation, with a depth similar to the real Herschel Science Demonstration Phase 
GOODS-N data taken as a part of HerMES. The palette is scaled from — 15mjy (white) to +20mjy (black). The panels show, from left 
to right, the raw, match-filtered, and PSF-filtered maps. The large circles indicate submm peaks with flux densities 5250 > lOmJy, and 
their radii are 3(Tr, where Cr is estimated as in Appendix |B 1 1 using a matching catalogue consisting of input sources with 5250 > 4mjy. 
The small circles and dots indicate input sources with flux densities 5250 > lOmJy, and 4mjy < 5250 < lOmJy, respectively. The 
match-filtered image is significantly less confused than the raw and PSF-filtered maps, but in many cases the submm peaks are still 
clearly blends of multiple sources. 



data, as expected given the relative beam sizes. Matched- 
filtering has an even greater potential to de-blend and im- 
prove the positional uncertainty for SPIRE than for BLAST 
due to the improved SNR. However, SPIRE peak catalogues 
will still generally represent blends of multiple sources. It 
will therefore be necessary to identify multiple potential 
counterparts to each peak, and we advocate performing a 
simultaneous fit to the data as we have done to the BLAST 
data in order to produce the least-biased photometry. Fi- 
nally, while not simulated here, the considerably larger cat- 
alogues of SPIRE peaks will enable estimates of priors for 
Likelihood Ratios (as in Sectionjsjl with considerably smaller 
statistical errors. In fact, LRs have already been use to pro- 



pose identifications for SPIRE sources in Smith et al. ( 2010 1 



However, they have calculated counterpart reliabilities fol- 



lowing Sutherland & Saunders ( 1992j ) , which we suspect may 
be biased due to the impact of confusion (see discussion in 
Appendix [b| . 



6 CONCLUSIONS 

We have performed a deep multi-wavelength study of indi- 
vidual peaks from the BLAST 250, 350, and 500 |J.m survey 
within the Extended Chandra Deep Field South (ECDF- 
S). By comparing the BLAST data with LABOCA ECDFS 
Submm Survey (LESS) maps at 870 \mi we have been able 
to greatly improve the positional uncertainties and longer- 
wavelength (Rayleigh- Jeans tail) SEDs for our sample. Com- 
pared to the earli er B LAST studies of|Dye et al.| ( |2009[ ), 
Ivison et al.l (120101) andlDunlop et al.|(|2010'), our methodol- 



ogy for identifying counterparts and measuring flux densities 
differ in several key respects: 

• It is recognized that peaks in the submm maps (par- 
ticularly the BLAST bands) are generally blends of several 
unrelated sources. It is therefore important to search for 
multiple counterparts to each peak in a matching catalogue. 



whereas earlier BLAST catalogue-based studies tended to 
focus on single counterparts. 

• We use a new 'matched filter' that compensates for 
source confusion when searching for submm sources. For 
the BLAST bands, in which the contribution to the total 
noise by source confusion is roughly a factor of two larger 
than the instrumental noise, there is a significant improve- 
ment in SNR of approximately 15-20%. In contrast, there is 
only a minimal improvement of about 5% for the LESS data 
since the confusion and instrumental noise components are 
approximately the same. 

• Identifications in our combined radio/IR matching cat- 
alogue have been made using Likelihood Ratios that incor- 
porate more prior information than traditional 'P' statistics. 
With a threshold spurious ID rate set at 10%, we find 52, 
50, and 31 matches to 64, 67, and 55 submm peak positions 
at 250, 350, and 500 [im, respectively. This is a significant 
improvement compared to the match rates obtained using 
10% cuts on P giving 45, 31, and 12 matches in the 24 jim 
catalogue, and 51, 35, and 17 matches in 1.4 GHz catalogue. 
Most of the gain is at 500 \im where the beam is so large that 
any additional information to help the process clearly makes 
a large difference. Combining the identifications made inde- 
pendently in each band, approximately 75% of the BLAST 
selected peaks have at least one potential counterpart. 

• We obtain submm photometry for the 118 unique iden- 
tifications in the matching catalogue by performing a simul- 
taneous fit of the submm PSFs at the precise ID locations in 
each of the four original (un-smoothed) submm maps. This 
procedure results in useful measurements for 76 sources; the 
remaining 42 sources are too close to one another to reliably 
disentangle the ffux densities produced by each of the con- 
tributing objects. In order to make further progress in such 
cases, prior information for the SEDs will likely be needed 
(see, for example, the work of Roseboom et al.|2009 l. 

• Our procedure for identifying counterparts and mea- 
suring their flux densities has compensated for much of the 
'flux-boosting' present in the original flux-limited submm 
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peak lists. Comparing the number counts of our final cat- 
alogue in the BLAST bands with the P{D) counts from 



Patanchon et al. ( 2009 1 we find good agreement above 
40 mjy at 250 [im, 30 mjy at 350 |J.m and with somewhat 
lower completeness reaching a maximum of ~ 66% above 
20 mJy at 500 ^m. 

We then proceeded to identify redshifts for the coun- 
terparts, and hence probe the rest-frame properties of our 
sample: 

• Of the 76 sources with usable submm photometry, 69 
counterparts have redshifts: 23 are optical spectroscopic red- 
shifts; 35 are optical photometric redshifts; and 11 are IRAC 
photometric redshifts. The median of the entire distribu- 
tion is 2: = 1.1 with an interquartile range 0.3-1.9. Restrict- 
ing ourselves to sub-samples above the flux-density limits 
mentioned previously shows a clear trend in increasing red- 
shift with observed wavelength: a median z = 1.1 with an 
interquartile range 0.2-1.9 at 250 ^.m; a median z = 1.3 
with an interquartile range 0.6-2.1 at 350 ii.m; and a median 
z = 1.6 with an interquartile range 1.3-2.3 at 500 (xm. In 
general, there is a higher-redshift tail for our sample than 



in the earlier study of Dye et al. ( 2009 1 which we believe is 



due to the fainter flux densities probed. 

• Fitting modified blackbody SEDs of the form Si, oc 
v^'^Bv{T) to the 69 sources with useable submm photome- 
try and redshifts, we establish rest-frame cold-dust temper- 
atures. The total distribution has a median T = 29K and 
interquartile range 23-36 K. There is almost no variation in 
the temperatures of the different fiux-limited sub-samples: 
a median T = 30 K with an interquartile range 25-39 K at 
250 [i.m; a median T = 30 K with an interquartile range 23- 
35 K at 350 |i.m; and a median T = 30 K with an interquartile 
range 22-38 K at 500 \i.m. These temperatures are systemat- 
ically warmer by about 7K than those inferred in [Dye et al.| 
( 2009 1 . While this discrepancy may be partly caused by the 



relative depths of the samples, we believe it is mostly due 
to improved compensation for confusion (which reduces the 
estimated flux densities at longer BLAST wavelengths more 
than at shorter wavelengths), and the inclusion of higher- 
resolution 870 vim data from LESS. 

• The primary improvement over earlier catalogue-based 
BLAST studies is the characterization of submm SEDs for 
the fainter, but higher-redshift [z ^ 1) and more-luminous 
galaxies. Even though the earlier analyses typically consid- 
ered single counterparts to the BLAST peaks, these coun- 
terparts are probably significant submm emitters (although 
the measured flux densities were typically high, especially 
in the most confused 500 vim band). Therefore the measured 
properties of the optical counterparts (approximate redshift 
distributions, galaxy types etc.), and even the FIR tempera- 
tures and luminosities of the brighter, lower-redshift galaxies 
are generally correct. 

• We have also fit star-forming galaxy SED templates 
from the library of [Dale et al.| ( |2001| | (which span the radio- 
near-IR) to the submm and FIR photometry of our sample. 
Generally speaking there is good correspondence between 
these models and the photometry in the radio and near/mid- 
IR, confirming that most of the emission is probably powered 
by star-formation. However, about 8% of the sample exhibits 
a significant excess in either (or both) of these wavelength 
regimes, which could indicate AGN activity. We note that 



these features are restricted to sources with luminosities pri- 
marily in the range 11.5 < log(LpiH,) < 12.8, near the bright 
end of our sample. 

• We compare the distribution of luminosity and cold 
dust temperature in our sample with the local-Universe mea- 



surement of Chapin et al. (20091. While the submm colours 



of our sources appear systematically cooler at all luminosi- 
ties compared to the local distribution, we have determined 
that the distribution is consistent with selection effects in 
our survey. We therefore find no evidence for evolution in 
the temperature-luminosity correlation out to z ~ 2.5. 

• Finally, we investigate the utility of the methods de- 
scribed in this paper for the new generation of Her- 
sc/ie//SPIRE surveys. Using a simple simulation we show 
that these new surveys are significantly confused, despite a 
factor of ~2 improvement in angular resolution over BLAST. 
We find that the matched-filter will yield a greater improve- 
ment in the detection of point sources than for BLAST due 
to the significantly lower instrumental noise of SPIRE. How- 
ever, even with these improvements, peaks in SPIRE maps 
will be confused, and we advocate the use of Likelihood Ra- 
tios to identify the expected multiple counterparts. 
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APPENDIX A: MATCHED FILTER 

Here we describe a filter for identifying point sources in 
maps containing significant contributions of both instrumen- 
tal noise and confusion due to blending of other point sources 
(see e.g. Tegmark & de 01iveira-Costa|1998 Barreiro et al 



[2003i |Vio et al.||2004[ for related studies of point sources m 
maps of the Cosmic Microwave Background). We formulate 



the problem as follows: we wish to find the 'matched filter', 
F, that maximizes the SNR one would obtain when cross- 
correlating F with the signal of interest, 5 - in our case, a 
point source whose shape is identical to the PSF - in the 
presence of noise, A''. We proceed by expressing the total 
SNR resulting from this operation, 0, in Fourier space, 



Signal 
Noise 



1/2 ■ 



(Al) 



In this expression the carets denote discrete Fourier trans- 
forms, and the index k runs over all elements in the spatial 
frequency domain. The superscript 'T' indicates that we use 
the transpose of F in the expression]^ Taking the partial 
derivatives of (f> with respect to each mode j in the filter, 
FJ, and setting them to 0, one can solve for the filter that 
maximizes the SNR: 



dFj 



(A2) 



1/2 



Fj\N,\' 



F 



T.,F^S, 



3/2 



(A3) 



In other words, in Fourier space the matched filter is sim- 
ply the PSF weighted by the inverse noise variance at each 
spatial frequency. In the case of an isolated point source in 
a field of white noise, this expression results in the PSF. 



since N is constant. This is the well-known result (e.g. Stet- 
son|1987 1 that the best way to find point sources in a noisy 
map is to convolve with the PSF. However, this procedure is 
only optimal if the noise is white. For our submm maps we 
consider two components of noise: instrumental white noise, 
A'w, and the confusion noise caused by other point sources. 

For the case at hand we estimate A'w (constant at all 
spatial frequencies) assuming typical noise values at the cen- 
tres of each submm ECDF-S map: 32, 29 , 27, and 2mJ>p]at 



250, 350, 500, and 870 ^m, respectively ( [Devlin et al.||2009 
Weiss et al.||2009 l. For the confusion noise, A'c, we use the 



counts inferred using P{D) analyses from [Patanchon et al.| 



( 2009 1 and Weiss et al. (20091, and knowledge of the PSF 



In the absence of noise (and assuming no spatial clustering 
of sources), the sky can be thought of as a random superpo- 
sition of PSF shapes with different amplitudes (i.e., 5 func- 
tions smoothed by the PSF and scaled by the point source 
flux densities). We therefore approximate the Fourier trans- 
form of the map using a white noise distributiorj^multiplied 

^ The cross-correlation of the filter with the signal is equivalent to 

convolving the signal with the complex conjugate transpose of the 

filter. Since we are dealing with real-valued signals this operations 

reduces to a transpose. The convolution is then a simple product 

in Fourier space between the two transformed quantities. 

^ These are the noise values in the raw un-smoothed maps at each 

wavelength. 

^ In practice, even if the angular locations of galaxies are un- 
clustered, the 5 function map power spectrum is not perfectly 
white because the flux densities of galaxies are drawn from the 
number counts rather than a Gaussian distribution. 
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by the Fourier transform of the transpose of the PSF (by the 
convolution theorem) ; in other words the power spectrum of 
confusion noise rolls off in the same way as the PSF, but 
with a different normalization. To determine how the power 
spectrum of the PSF should be scaled we simulate maps in 
each band as described above, and use the ratios of the stan- 
dard deviations of these noise realizations to the standard 
deviations of the PSFs. 

As an example of this procedure, the estimated noise 
power spectrum at 500 (xm, and resulting matched filter in 



real space (both azimuthally averaged) are shown in Fig. Al 



One sees that the matched filter is somewhat narrower than 
the PSF filter (which has been used in most previous SMG 
studies) and has ringing at larger angular scales, giving it 
similar properties to the 'Mexican hat' kernel (see for exam- 



ple discussion and references within Barnard et al. 20041 



In this particular case, the first large negative ring effec- 
tively removes a local baseline estimated at radii ~ 50 arcsec, 
and hence corrects the flux density slightly for other nearby 
blended sources. The matched filter is also qualitatively sim- 
ilar to the filter used to detect sources in deep AzTEC maps 
(e.g., Scott et al.|2008 Perera et al.^2008J . However, in those 
cases only sources of atmospheric and instrumental noise 
are considered, since the spatial noise power spectrum is 
estimated using 'jackknife maps' (produced by differencing 
alternating portions of the data), which explicitly removes 
the contribution of astronomical sources. We emphasize that 
here we are using such a filter to optimally extract sources 
in a confused background. Note also that in a language fa- 
miliar in other fields, our filter is essentially a Wiener filter 
which includes source confusion explicitly as a noise term. 

In addition to the calculation we describe here, we also 
convolved simulated maps (containing both point sources 
and instrumental noise) with Gaussian filters of different 
widths, and determined the FWHMs that optimized the 
SNR of point sources. These simple tests gave results com- 
parable to the more detailed calculation above, although the 
SNR was slightly lower due to the lack of ringing to com- 
pensate for the noise at larger angular scales. 

It should be noted that in the development of our fil- 
ter we have assumed that the instrumental noise is constant 
across the map. Of course, in general this is not true, and 
one could in principle derive a different filter at each posi- 
tion. For the BLAST maps under discussion, a separate filter 
should certainly be used to detect sources in BGS-Wide and 
BGS-Deep, but for the present study the noise was close 
to uniform across the ECDF-S, and hence we used a single 
filter for the whole map. 

As a consequence of Eq. |A3| the peak SNR for an iso- 
lated point source in a field of white noise (for which Nk 
is constant) is obtained by convolution with a Gaussian of 
the same FWHM as the instrumental PSF. In the oppo- 
site extreme, for a map of confused point sources with no 
instrumental noise, the optimal filter is the inverse of the 
PSF in Fourier space, i.e., the map is de-convolved by the 
beam. For real data in which instrumental noise and con- 
fusion are both present, the degree to which the map is ac- 
tually de-convolved depends on the relative amplitudes of 



the two noise components (see left panel of Fig. All. The 
source confusion depends on the number counts (and there- 
fore the SEDs and redshift distribution) of galaxies in the 
observed band, and beam size, and the instrumental noise 



depends on the detector sensitivity and integration time. For 
the data described here, these two contributions are compa- 
rable, since the submm surveys of the ECDF-S were de- 
signed to be just confusion-limited. The bottom-line is that 
in the BLAST bands, considering both noise terms, con- 
volving with kernels that are smaller than the instrumen- 
tal PSFs can give approximately 15-20% improvement in 
the SNR — quite significant considering the majority of the 
sources under discussion have instrumental SNR of about 
4(T (see Table [l|! The improvement is less impressive in the 
LESS 870 \im image, since the instrumental noise dominates 
the RMS resulting from confusion. Most previous studies of 
SMGs in the regime where confusion is important have ex- 
tracted sources using the PSF, and our results show that 
they have smoothed away some of the point source informa- 
tion in the maps. 



APPENDIX B: LIKELIHOOD RATIOS 

Here we describe our formulation of Likelihood Ratios (LR) 
as a means for identifying counterparts to submm peaks. 
The matching catalogue for which we estimate priors is de- 
scribed in Section [2.2| To begin, we make two basic assump- 
tions: (i) that the radial offset of a potential ID is uncorre- 
lated with its other properties; and (ii) that the submm fiux 
density is uncorrelated with its other properties. The first 
assumption is fairly standard and uncontroversial, although 
the latter could be considered problematic. For example, the 
brightest peaks at 250 vim are known to be low-luminosity 
IRAS galaxies, whereas the fainter sources are thought to 
be a mixture of both low- and high-redshift galaxies. How- 
ever, for now we take the practical route and choose not to 
make any distinction based on these properties so that we 
can estimate priors with reasonable SNR given the data. We 
then express the differential density of true counterparts to 
submm peaks as a function of their offsets, r, flux densi- 
ties, S (both at 24 [im and 1.4 GHz), and colour c (IRAC 
[3.6] - [4.5]): 



nc{S,c,r)dS dcdr — q{S, c)f{r)dSdcdr, 



(Bl) 



where q{S,c) is the distribution of counterpart flux densities 
and colours (in the matching catalogue), and /(r) is the po- 
sitional probability distribution as a function of radial offset 
r. We assume a symmetric Gaussian probability distribution 
as a function of orthogonal positional coordinates which re- 
sults in the Rayleigh radial probability distribution. 



f{r) 



(B2) 



The normalization is chosen so that f{r)dr = 1. Follow- 
ing Sutherland & Saunders ( 1992[ ) q{S,c) is normalized so 
that it integrates to the average expected number of coun- 
terparts per submm peak. However, while those authors as- 
sumed that this number is in the range 0-1 (i.e., the emis- 
sion is produced by a single source, which may or may not 
be present in the matching catalogue), we consider cases in 
which the number of counterparts may be greater than this. 

We can also estimate the differential density of back- 
ground sources as a function of these same properties 

n^{S,c,r)dSdcdr = 27vrp{S, c)dS dcdr, (B3) 

where p{S, c) is the surface density of background sources as 
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Figure Al. left: The peak- normalized total azimuthally averaged 500 |xm angular noise power spectrum (solid line) and its two compo- 
nents: instrumental (white) noise (dotted line), and confusion from point sources (dashed line). The latter function has the same shape 
as the PSF, since point sources are modelled as a spatially unclustered collection of point sources smoothed by the PSF. The roll-off at 
K < 0.003 in the PSF is a result of the spatial whitening filter that has been applied to the BLAST maps pevlin et ar][2009j . right: 
The azimuthally averaged 500 |J.m matched filter (solid line) compared with the PSF (dashed line), now in angular rather than frequency 
space. The matched filter is obtained by dividing the power spectrum of the PSF (same shape as the dashed line in the left-hand plot) 
by the noise power spectrum (solid line in the left-hand plot), and taking the inverse Fourier Transform. For the noise sources considered 
here, this procedure augments high-frequency components of the PSF, resulting in a narrower central profile, while introducing ringing 
at angular scales ^ 50 arcsec. 



a function of S and c, and multiplying by 27rr converts this 
quantity to the infinitesimal number density of sources at a 
distance r from the submm peak. 

Therefore, given S, c and r for the jth candidate coun- 
terpart to the iih submm peak, the relative number of ex- 
pected true counterparts to background sources is the LR: 



Li,j 



rih 



(B4) 



2Tva?p(Sj,Cj) 

Candidates with large values of Lij are more hkely to be 
associated with the submm emission. 

Note that while the LR contains all of the information 
about the potential identification, it is not itself a probabil- 
ity distribution. In the case that only a single source in the 
matching catalogue is believed to correspond to the submm 
peak, and where q{S,c), p{S,c), and a are all known pre- 
cisely, it is possible to calculate the reliability R, the proba- 
bility that the candidate is the unique counterpart, following 



Sutherland & Saunders (19921. Using this technique on a 



case for which there are multiple strong candidate matches, 
for example, the reliability is divided among them in such a 
way that their sum does not exceed 1. For the case at hand, 
since we consider the possibility of multiple real counter- 
parts to a single submm peak, and since our estimates of 
the priors are ultimately quite noisy, Eq. 5 from [Sutherland] 
I&: Saunders| ( |1992[ ) does not apply. Instead we will take the 
approach of using Monte Carlo simulations to establish a 
threshold LR that will provide candidates with a desired 
false identification rate. 



Bl Background Counts and Positional 
Uncertainties 

We estimate priors directly from the data themselves. First, 
p{S,c), the background number counts, are determined by 



binning the entire matching catalogue (Section 2.2 1 as a 
function of the two flux densities and colour, and dividing 
by the survey area. This should give a good estimate for the 
average properties of the background sources since the en- 
tire catalogue contains 9216 entries, whereas there are only 
of order ~ 60 submm peaks in each band. 

The positional uncertainty parameter, a^, is more dif- 
ficult to estimate from the data. We proceed by counting 
the number of catalogue sources around submm peaks, and 
checking for statistically significant excesses compared to 
the background counts predicted by p{S,c), as a function 
of search radius. The expectation, if each submm peak were 
produced by a single source in the matching catalogue, is 
that the cumulative excess would grow from at a search 
radius of 0, and converge to 1 at large radii. 

We show the differential excess distribution (calculated 



in annuli of difi'erent sizes) for each submm band in Fig. Bl 



Uncertainties for each data point are estimated from Poisson 
counting statistics (considering both the number of objects 
in the full catalogue for the uncertainties in the background 
counts, and the number of objects in the measurement an- 
nulus). We fit the Rayleigh radial offset distribution to the 
excess counts, multiplying the single-parameter expression 
for /(r) in Eq. |B2| by a second parameter, E, giving the 
integrated excess number of sources around submm peaks 
(note that E = f f q{S,c)dSdc). The smooth models are 
shown as solid lines in the differential excess plots. In addi- 
tion, the solid histograms show the models integrated across 
the same bins as the data, giving the model predictions used 
to calculate x^- Clearly, E is significantly larger than 1 in 
each submm band, demonstrating that the submm peaks 
from our survey are typically blends of several sources in 
the matching catalogue. One could use this approach to de- 
termine a criterion which might be called 'counterpart con- 
fusion' - for a given source map and catalogue (in a different 
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waveband) one could conclude that the source identification 
process will be relatively straightforward if £ < 2 (say) , but 
significantly complicated by confusion ii E ^ 2. 



B2 Flux density and colour priors 

Next we estimate the joint distribution of flux density and 
colour for counterparts, q{S,c). Since we assume that these 
properties are un-correlated with their distance from the 
submm centroid, we first identify the search radius within 
which the SNR of the excess counts is highest, r^. As this 
radius is increased from 0, the background counts grow as 
Tg , but the counts around submm positions grow even faster 
as the counterparts are included. The ratio of the difference 
(the excess), compared to the Poisson uncertainty in the dif- 
ference, therefore increases until the number of new coun- 
terparts drops significantly. We find that for our matching 
catalogue and submm peaks lists the maximum excess SNR 
are achieved at 23, 23, 35 and ISarcsec for 250, 350, 500 and 
870 [im, respectively. Using these search radii, we then com- 
pare the normalized histogram of 24 jim and 1.4 GHz flux 
densities for the entire matching catalogues with the nor- 
malized histogram of excess sources around submm peaks, 



p{S24) and p{Si), in Figs. B2 and B3 This operation shows 
that the extra sources around submm peaks tend to be 
brighter at 24 (xm and 1.4 GHz on average than the entire 
populations in the BLAST bands. 

In a similar way we compare the log j^q (5*3. 6/54. 5) colour 
of the background population to the excess around submm 
positions. This IRAC colour was chosen because it is the pri- 



mary discriminator for redshift used in Devlin et al. (20091, 
[Marsden et al.| ([2009'), and [Pascale et al.| ( |2009[ ) - with 
the trend that smaller ratios (i.e. redder colours) correlate 
with higher redshifts. Earlier studies have also used IRAC 
colours as a crude redshift estimator (e.g., Pope et al.|2006 



Yun et al.|2008| [Wilson et al.l2008l ). The normalized excess 
colour distribution, p(c) (Fig. |B4[ ) agrees with this expecta- 
tion. There is a trend from bluer colours starting at 250 vim 
(submm excess and background colours nearly indistinguish- 
able), to significantly redder colours at 500 and 870 micron. 

Individually, p{S24), p{St), and p{c) give normalized es- 
timates of g(524, Sy, c) marginalized over the remaining vari- 
ables. If they were completely independent of one another, 
we could estimate q{S, c) as 



q{S,c)~Exp{S24)p{Sr)pic), 



(B5) 



where E is the total excess measured in Fig. |B1[ Ideally 
we would like to bin the excess counts in cubes of 5*24, Si, 
and c simultaneously. Unfortunately, given the sample size, 
it is not possible to obtain a statistically significant mea- 
surement of this distribution. Similarly, we approximate the 
background counts 



p{S,c)^p{S24)p{Si)pic) 



(B6) 



We proceed under the assumption that these quantities are 
independent, but in the next section use Monte Carlo simu- 
lations to establish a reasonable threshold for this approxi- 
mate LR to obtain counterparts. 



B3 Normalized LR 

Since each submm peak appears to be produced by a blend 
of several matching catalogue sources, we cannot use the 



Sutherland & Saunders 1 1992 1 reliability, R, as a normalized 



probability that a given candidate is the single counterpart. 
Furthermore, its absolute normalization depends on precise 
measurements of a^, p{S,c) and q{S,c), all of which have 
moderate uncertainties for this sample (especially q{S, c), as 
noted in the previous section). 

Instead, we compare the LRs for all potential candidate 
matches to submm peaks, with the distribution of LRs for 
false matches to 10,000 random positions. This procedure 
enables us to select a threshold LR that gives an acceptable 
rate of false-positives, which we set to 10%. This approach 
to identifying the threshold LR is similar to that of |Mann| 
et al.| (19971. The results of our calculation are shown in 



Fig. |B5| The solid histogram is the distribution of LR for all 
matching catalogue sources out to 60 arcsec (this maximum 
search radius easily contains > 99% of the true counterparts) 
from each submm position. The dotted histograms show the 
'background' distribution obtained by searching for counter- 
parts around random positions. As expected, for large values 
of the LR, there are more sources around submm peaks than 
random positions. Also, unsurprisingly, the distinction be- 
tween the two curves is greatest at 870 \im (which has the 
smallest positional uncertainties), and progressively worse 
across the BLAST bands to longer wavelengths. 

A threshold in the LR is chosen such that the fraction 
of sources with larger LRs for the simulation (random po- 
sitions) is 10%. This cut is indicated as a vertical dashed 
line in Fig. |B5[ For convenience, we normalize the LRs by 
these threshold values, so that potential counterparts have 
values that are greater than 1. Also, noting the ratio of the 
density of sources around submm and random positions, a 
candidate counterpart at this threshold is about three times 
more likely to be real than spurious (this value does not, 
however, reflect on the absolute probability that the source 
is either real or spurious). 

The results of applying this cut to potential identifica- 
tions around submm positions are given in Table |B1| For 
each band the total number of expected counterparts is the 
product of the number of submm positions with the average 
excess of matching catalogue sources around each spot. The 
estimated number of true IDs found is then the number of 
sources detected above the cut on the LR, after subtract- 
ing the 10% spurious fraction. This test shows that as many 
as 32% of the individual matching catalogue sources con- 
tributing to the 870 vim peaks have been identified, and as 
few as 13% at 500 vim. The trend of this result is certainly 
expected, due to the increased confusion caused by larger 
beams. However, the extent of the problem of identifying 
confused sources is larger than some would have expected. 



APPENDIX C: DATA TABLES 

The following data tables are provided in this Appendix: 
the match-filtered submm peak lists in the BLAST bands 
(Tables C1|C3 1, as well as the LESS peaks from Weiss et al 



(2009) that land within the region that was analyzed (Ta- 



ble C4 1 ; the correspondence between the matching catalogue 
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Figure Bl. Excess source counts in annuli around submm peaks compared to the expected background as a function of search radius. 
Poisson uncertainties are plotted. The soHd Hnes are fits of a Gaussian radial probability density function p(r) = (Er/a^) exp(— r^/2o-^), 
where the scale factor E is the total excess counts encountered on average (fit values are indicated in each panel). This function is fit to 
the differential excess counts (since the uncertainties are uncorrelated) . The histograms in the differential plots show the smooth model 
integrated across each bin (giving the actual predicted model values used to calculate X?)- 



Table Bl. Excess count and ID statistics. The number of peaks are from the BLAST and LESS catalogues, using 3.75(j and 3.72a 
thresholds, respectively. The excess counts per source and ffr are the best-fit positional uncertainty model parameters (with Xr giving 
the reduced chi-squared in each case) — see Section |B1[ The fact that the excesses are greater than 1 is an indication of how confused 
the data are, and a warning that finding single counterparts is challenging. Expected numbers of IDs are calculated from the product 
of the number of submm peaks with the average excesses per source. 'Found IDs' are the total number of potential counterparts with 
false-identification rates of 10% in the matching catalogue. The 'Real ID' rate subtracts the expected number of false IDs for the list 
(10% of the number of submm peaks; see discussion of threshold LR values in Section|3]l, and also expresses that quantity as a percentage 
of the total expected. 
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t There are 81 peaks from LESS that land in the region of coverage considered here, but only 42 peaks are cross-matched to the BLAST 
peak lists. The full list of 81 peaks is used to estimate priors even though we only discuss the properties of the smaller subset of 

cross-matched peaks in this paper. 
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Figure B2. MIPS 24 \xm flux density distributions of excess sources around submm pealts in all four bands (solid histogram with 
Poisson uncertainties), our estimate of p{S24,), compared to the general population (dotted histogram), our estimate of p{S24,)- The bin 
at contains all radio sources that are not present in the FIDEL 24 |^m catalogue. Galaxies selected near submm positions are on average 
brighter at 24 (xm than galaxies selected randomly, and this is true for all four bands. 



and each of the submm lists (Tables C5 1; and the properties 
of the matched sources (re-measured submm flux densities, 
SED fits and redshifts. Table [C6 1. 



APPENDIX D: POSTAGE STAMPS 



This Appendix contains submm, IRAC and radio postage 
stamps for all 118 proposed identifications to the submm 
peaks (Fig. Dl |. 
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Figure B3. VLA 1.4 GHz flux density distributions of excess sources around submm peaks in all four bands (solid histogram with 
Poisson uncertainties), our estimate of p(S'i,4), compared to the general population (dotted histogram), our estimate of p(Si.4). The bin 
at contains all 24 i^m FIDEL sources that do not exhibit significant radio flux. Galaxies selected near submm positions are on average 
brighter at 1.4 GHz than galaxies selected randomly, and this is true for all four bands. We also note that many more of the radio sources 
with 24 |xm emission are found near submm positions than random radio sources. Many of these radio sources without 24 ]x.ra emission 
could be spurious given the low significance cut on the catalogue that has been used (see Section |2.2.2| l. Our prior estimation procedure 
will give much greater weight to potential radio identifications that also exhibit 24 yira emission. 
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Figure B4. IRAC logj(,(53.6/S'4.5) colour distributions of excess matching catalogue sources around submm positions in all four bands 
(solid histogram with Poisson uncertainties), our estimate of p(c), compared to the general population (dotted histogram), our estimate 
of p(c). There is a clear trend of redder colours in the IRAC bands for galaxies selected at longer submm wavelengths. 
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Figure B5. Distribution of likelihood ratios (LR) for potential counterparts in the matching catalogue out to 60arcscc from submm 
peak positions (solid histograms) compared to random positions (dotted histograms). Thresholds for identifying potential counterparts 
are given by the vertical dashed lines, such that the integrated tail of the distribution for random positions gives a false positive rate of 
10%. All LRs are normalized to this value (i.e., potential IDs with LR > 1.00 have less than a 10% chance of being false). 
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Table CI. 250 ^im peak list produced from matched-filtered maps that land within the coverage of the matching catalogue (Figure [2|. 
On average, each of these peaks is a blend of 3.2 sources from our matching catalogue. The 'Previously Published Name' is taken from 
the supplement to [Devlin et al.| l |2009| l when available (new peaks have no entry in this column), but the short submm 'SID' is used 
throughout this paper. The 'Match' column refers to IDs in Table [C5] The positions are centroids of the submm peaks. The flux densities 
are raw values from the maps (with no correction for flux boosting), and the noises and SNR refer to instrumental noise only (excluding 
source confusion). 
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91 


BLAST 


J033223- 


273642 


(250n.m) 


59 


03 


32 


22 


81 


-27 


36 


43 


4 


66 


± 


14 


4.72 


92 


BLAST 


J033140- 


275633 


(250n.m) 


18 


03 


31 


40 


33 


-27 


56 


36 


3 


66 


± 


14 


4.70 


93 


BLAST 


J033211- 


275859 


(250n.m) 




03 


32 


11 


86 


-27 


59 


00 


8 


66 


± 


14 


4.70 


97 


BLAST 


J033222- 


280019 


(250n.m) 


58 


03 


32 


22 


14 


-28 


00 


21 


1 


65 


± 


13 


4.66 


99 


BLAST 


J033129- 


275910 


(250n.m) 


9 


03 


31 


29 


19 


-27 


59 


11 


1 


66 


± 


14 


4.64 


111 


BLAS 1 


J033130- 


275604 


(250n.m) 


7,8 


03 


31 


29 


99 


— 27 


56 


02 


1 


63 


± 


14 


4.53 


112 










63 


03 


32 


29 


95 


-27 


43 


14 


4 


64 


± 


14 


4.52 


113 


BLAST 


J033235- 


274932 


(250n.m) 


67 


03 


32 


35 


24 


-27 


49 


28 


2 


62 


± 


13 


4.51 


121 


BLAST 


J033230- 


275905 


(250n.m) 


64 


03 


32 


30 


17 


-27 


59 


05 


9 


63 


± 


14 


4.44 


130 


BLAST 


J033259- 


273535 


(250n.m) 


93 


03 


32 


59 


71 


-27 


35 


33 


4 


61 


± 


14 


4.37 


146 


BLAST 


J033151- 


274431 


(250n.m) 


31,32 


03 


31 


51 


13 


-27 


44 


38 


1 


60 


± 


13 


4.29 


151 


BLAST 


J033147- 


274147 


(250n.m) 




03 


31 


47 


89 


-27 


41 


41 


2 


60 


± 


14 


4.27 


154 












03 


32 


06 


88 


-27 


38 


00 


7 


60 


± 


14 


4.25 


178 


BLAST 


J033232- 


275304 


(250n.m) 




03 


32 


33 


03 


-27 


53 


04 


8 


59 


± 


14 


4.18 


185 


BLAST 


J033317- 


274118 


(250n.m) 


110 


03 


33 


17 


38 


-27 


41 


12 


5 


59 


± 


14 


4.15 


188 


BLAST 


J033154- 


274406 


(250n.m) 


37 


03 


31 


54 


28 


-27 


44 


04 


2 


58 


± 


14 


4.15 


210 


BLAST 


J033135- 


274705 


(250n.m) 




03 


31 


35 


03 


-27 


47 


09 


1 


57 


± 


14 


4.09 


222 


BLAST 


J033200- 


280234 


(250n.m) 




03 


32 


00 


67 


-28 


02 


32 


1 


57 


± 


14 


4.06 


226 












03 


31 


29 


94 


-27 


39 


41 


9 


58 


± 


14 


4.05 


235 


BLAST 


J033156- 


280306 


(250n.m) 




03 


31 


56 


85 


-28 


03 


09 


4 


57 


± 


14 


4.03 


244 










17 


03 


31 


39 


58 


-27 


41 


29 


5 


56 


± 


14 


4.00 


246 


BLAST 


J033211- 


280242 


(250n.m) 


45 


03 


32 


11 


14 


-28 


02 


42 


6 


56 


± 


14 


4.00 


250 


BLAST 


J033251- 


274530 


(250n.m) 




03 


32 


51 


55 


-27 


45 


33 


3 


56 


± 


14 


3.99 


255 


BLAST 


J033144- 


275521 


(250n.m) 


24 


03 


31 


44 


76 


-27 


55 


19 


5 


55 


± 


13 


3.98 



BLAST and LABOCA observations of the ECDF-S 



SID Previously Published Match R.A. Dec. Flux Density SNR 



250 




Name 




ID 


(■^ 


m 




(° 




") 


(mjy) 






279 


TIT \ orn 


J033311- 


-274313 


(250 |xm) 




03 


33 


11, 


.13 


—27 


43 


12.0 


56 


± 


14 


3. 


.93 


289 


TIT A orn 

BLAbi 


J033217- 


-275054 


(250|J.m) 


54 


03 


32 


17, 


.43 


—27 


50 


55.7 


54 


± 


13 


3. 


.91 


298 


BLAST 


J033243- 


-275146 


(250|xm) 


79,80 


03 


32 


43, 


.50 


-27 


51 


45.6 


54 


± 


14 


3, 


,89 


312 


BLAST 


J033228- 


-273545 


(250|J.m) 




03 


32 


28, 


.66 


-27 


35 


53.0 


54 


± 


14 


3, 


,86 


316 


BLAST 


J033238- 


-275651 


(250(J.m) 


72 


03 


32 


37, 


.99 


-27 


56 


57.6 


53 


± 


13 


3, 


,86 


325 


BLAST 


J033225- 


-273822 


(250|xm) 


60 


03 


32 


24, 


.97 


-27 


38 


22.2 


54 


± 


14 


3, 


,84 


328 


BLAST 


J033205- 


-280054 


(250|J.m) 


42 


03 


32 


05 


.54 


-28 


00 


54.2 


53 


± 


13 


3, 


,84 


348 


BLAST 


J033150- 


-275100 


(2501^111) 




03 


31 


50 


.24 


-27 


51 


11.0 


53 


± 


13 


3, 


,81 


354 


BLAST 


J033248- 


-274443 


(250^m) 




03 


32 


47, 


.55 


-27 


44 


58.3 


53 


± 


14 


3. 


.80 


358 


BLAST 


J033210- 


-274253 


(250^m) 


47 


03 


32 


10, 


.51 


-27 


42 


55.3 


53 


± 


14 


3. 


.80 


369 










87 


03 


32 


49, 


.09 


-27 


36 


19.7 


53 


± 


14 


3. 


.78 


378 


BLAST 


J033131- 


-274601 


(250Hm) 




03 


31 


31 


.77 


-27 


46 


05.1 


52 


± 


14 


3, 


,76 


381 


BLAST 


.1033311- 


-275226 


(250|J.m) 




03 


33 


11 


.65 


-27 


52 


27.0 


52 


± 


13 


3, 


,76 


389 


BLAST 


J033251- 


-275936 




88 


03 


32 


51, 


.85 


-27 


59 


40.9 


53 


± 


14 


3, 


,75 
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Table C2. 350 |^m peak list produced from matched-filtered maps. Columns have the same meaning as in Table [CT] On average, each 
of these peaks is a blend of 3.4 sources from our matching catalogue. 



SID Previously Published Match R.A. Dec. Flux Density SNR 



350 




Name 




ID 




m 




(° 




") 


(mjy) 






3 


BLAST 


J033234- 


275531 


(350n.m) 


66 


03 


32 


34 


88 


-27 


55 


31 


3 


89 




10 


8 


24 


7 


BLAST 


J033249- 


275833 


(350n.m) 


86 


03 


32 


49 


49 


—27 


58 


33 


8 


79 


± 


10 


7 


36 


8 


BLAST 


J033229- 


274414 


(350n.m) 


61,62 


03 


32 


29 


37 


-27 


44 


16 


4 


74 




10 


6 


95 


10 


BLAST 


J033220- 


275647 


(350n.m) 




03 


32 


20 


15 


—27 


56 


45 


4 


73 




10 


6 


82 


11 


BLAST 


J033258- 


274328 


(350n.m) 


92 


03 


32 


58 


39 


-27 


43 


25 





72 




10 


6 


73 


12 


T AO rri 

BLAST 


J033229- 


274305 


(350n.m) 


63 


03 


32 


29 


49 


—27 


43 


03 


4 


71 




10 


6 


72 


14 


BLAST 


J033328- 


275700 


(350n.m) 


118 


03 


33 


27 


65 


—27 


57 


03 


4 


72 


± 


10 


6 


66 


16 


T AO rn 

BLAST 


J033129- 


275722 


(350n.m) 


10,11 


03 


31 


29 


61 


—27 


57 


22 


6 


70 




11 


6 


41 


19 


BLAST 


J033204- 


274650 


(350n.m) 


43 


03 


32 


04 


59 


-27 


46 


51 


7 


67 




10 


6 


35 


20 


BLAST 


J033207- 


275815 


(350n.m) 


44 


03 


32 


07 


34 


-27 


58 


18 


1 


68 




10 


6 


30 


22 


BLAST 


J033128- 


273927 


(350n.m) 


2,3,4 


03 


31 


28 


02 


—27 


39 


32 


2 


71 




11 


6 


25 


31 


BLAST 


J033311- 


274135 


(350n.m) 


105 


03 


33 


11 


16 


—27 


41 


32 


5 


64 


± 


10 


5 


95 


34 


BLAST 


J033247- 


274224 


(350n.m) 


83 


03 


32 


47 


24 


-27 


42 


31 


9 


62 




10 


5 


86 


35 


BLAST 


J033150- 


274333 


(350n.m) 


30 


03 


31 


49 


77 


-27 


43 


28 


2 


63 




10 


5 


84 


37 


BLAST 


J033140- 


274435 


(350n.m) 


19,20 


03 


31 


40 


71 


—27 


44 


36 


9 


62 


± 


10 


5 


76 


42 


BLAST 


J033138- 


274122 


(350n.m) 


17 


03 


31 


39 


04 


—27 


41 


22 


3 


60 


± 


10 


5 


67 


43 


BLAST 


J033318- 


274606 


(350n.m) 


111,112 


03 


33 


18 


39 


-27 


46 


13 


4 


60 


± 


10 


5 


64 


45 


BLAST 


J033210- 


275206 


(350n.m) 


46 


03 


32 


11 


53 


-27 


52 


03 


8 


61 


± 


10 


5 


62 


47 


BLAST 


J033151- 


274428 


(350n.m) 


31,32 


03 


31 


51 


20 


—27 


44 


33 


4 


59 


± 


10 


5 


56 


49 


BLAST 


J033135- 


275448 


(350n.m) 


13,14,15 


03 


31 


35 


44 


—27 


54 


47 


5 


60 


± 


10 


5 


56 


52 


T AO rn 

BLAST 


J033252- 


273756 


(350n.m) 




03 


32 


52 


12 


—27 


37 


47 


9 


58 




10 


5 


36 


58 


BLAST 


J033321- 


275513 


(350n.m) 


117 


03 


33 


21 


74 


-27 


55 


13 


1 


56 




10 


5 


19 


64 


BLAST 


J033302- 


275640 


(350n.m) 


95 


03 


33 


02 


20 


—27 


56 


43 


7 


54 


± 


10 


5 


05 


65 


BLAST 


J033218- 


275207 


(350n.m) 


55,56 


03 


32 


18 


12 


—27 


52 


12 


5 


54 


± 


10 


5 


03 


69 


BLAST 


J033237- 


273541 


(350n.m) 


70,71 


03 


32 


37 


62 


—27 


35 


42 


4 


54 




11 


4 


96 


79 


BLAST 


J033307- 


275833 


(350n.m) 




03 


33 


06 


98 


-27 


58 


35 


5 


52 




10 


4 


84 


80 


BLAST 


J033217- 


275906 


(350n.m) 


53 


03 


32 


17 


78 


-27 


59 


20 


2 


52 


± 


10 


4 


83 


84 


BLAST 


J033318- 


274131 


(350n.m) 




03 


33 


19 


72 


-27 


41 


38 


7 


52 


± 


10 


4 


79 


85 


BLAST 


J033230- 


275905 


(350n.m) 


64 


03 


32 


30 


00 


-27 


59 


02 


4 


51 




10 


4 


79 


88 


BLAST 


J033211- 


273731 


(350n.m) 


49 


03 


32 


11 


34 


-27 


37 


24 


1 


51 




10 


4 


76 


90 


BLAST 


J033253- 


280107 


(350n.m) 


89 


03 


32 


53 


05 


—28 


01 


11 


3 


51 




10 


4 


76 


91 


BLAST 


J033321- 


273908 


(350n.m) 




03 


33 


21 


00 


—27 


39 


10 





51 




10 


4 


75 


92 


BLAST 


J033308- 


275130 


(350n.m) 


104 


03 


33 


08 


61 


-27 


51 


29 


3 


50 




10 


4 


74 


93 


BLAST 


J033153- 


274952 


(350n.m) 


36 


03 


31 


54 


84 


-27 


49 


37 


7 


51 




10 


4 


73 


95 


BLAST 


J033152- 


280325 


(350n.m) 


33 


03 


31 


53 


14 


-28 


03 


35 


7 


51 




10 


4 


72 


98 










108 


03 


33 


15 


41 


—27 


45 


32 


9 


50 




10 


4 


69 


108 


BLAST 


J033215- 


273929 


(350n.m) 


51 


03 


32 


15 


25 


-27 


39 


31 


8 


50 




10 


4 


62 


112 


BLAST 


J033316- 


275056 


(350n.m) 


109 


03 


33 


17 


16 


-27 


50 


55 


7 


48 


± 


10 


4 


60 


115 


BLAST 


J033305- 


274414 


(350iJ.m) 


96 


03 


33 


06 


26 


-27 


44 


14 


1 


49 




10 


4 


59 


116 


BLAST 


J033245- 


275739 


(350n.m) 


81 


03 


32 


46 


04 


-27 


57 


35 


2 


49 




10 


4 


59 


117 


BLAST 


J033134- 


274629 


(350n.m) 




03 


31 


34 


56 


-27 


46 


25 


9 


49 


± 


10 


4 


57 


124 


BLAST 


J033153- 


273931 


(350n.m) 


34,35 


03 


31 


52 


62 


-27 


39 


29 


9 


49 




10 


4 


54 


128 


BLAST 


J033243- 


275509 


(350n.m) 


76,77 


03 


32 


42 


86 


-27 


55 


13 


1 


48 




10 


4 


51 


137 


BLAST 


J033217- 


275906 


(350n.m) 


52 


03 


32 


17 


89 


-27 


58 


46 


5 


47 




10 


4 


44 


144 


BLAST 


J033152- 


275035 


(350n.m) 




03 


31 


51 


22 


-27 


50 


41 


8 


46 




10 


4 


40 


145 


BLAST 


J033308- 


274813 


(350n.m) 


100,101,103 


03 


33 


08 


49 


-27 


48 


16 


5 


47 




10 


4 


40 


152 


BLAST 


J033145- 


274153 


(350n.m) 




03 


31 


44 


98 


-27 


41 


50 


5 


46 


± 


10 


4 


32 


154 


BLAST 


J033238- 


274620 


(350n.m) 




03 


32 


38 


52 


-27 


46 


18 


5 


46 




10 


4 


32 


161 


BLAST 


J033212- 


275459 


(350n.m) 




03 


32 


12 


74 


-27 


54 


59 


1 


45 




10 


4 


29 


169 


BLAST 


J033225- 


273810 


(350n.m) 


60 


03 


32 


24 


93 


-27 


38 


17 


3 


46 




10 


4 


25 



BLAST and LABOCA observations of the ECDF-S 



SID Previously Published Match R.A. Dec. Flux Density SNR 



350 




Name 




ID 


0" 


m 


^) 


(° 




") 


(mjy) 






170 


TIT A orri 


J033201- 


-274142 


(350|xm) 


41 


03 


32 


02 


.02 


—27 


41 


38.0 


45 


± 


10 


4. 


,25 


177 


OT A OT^ 

oLAb i 


J033129- 


-275911 


(3o0|xmj 


9 


03 


31 


29, 


.63 


—27 


59 


09.1 


46 


± 


10 


4. 


.22 


184 












03 


32 


10 


.66 


—27 


58 


56.9 


45 


± 


10 


4. 


,21 


187 


T) T A cn 

BLAb 1 


J033141- 


-275529 


{350(J.m) 


23 


03 


31 


41, 


.76 


— 27 


55 


28.4 


44 


± 


10 


4, 


18 


188 


T) T AC 

BLAb 1 


J033147- 


-274755 


(350(J.m) 


28 


03 


31 


47, 


.60 


— 27 


47 


52.9 


45 


± 


10 


4, 


,18 


199 


BLAST 


J033217- 


-275407 


(350(J.m) 




03 


32 


17, 


.92 


-27 


54 


08.0 


44 


± 


10 


4, 


,13 


235 


BLAST 


J033127- 


-274430 


(350(J.m) 


1 


03 


31 


27, 


.90 


-27 


44 


32.8 


43 


± 


10 


4, 


,04 


259 


BLAST 


J033235- 


-280135 


(350|J.m) 


68,69 


03 


32 


35, 


.90 


-28 


01 


48.7 


42 


± 


10 


3, 


,97 


264 


BLAST 


J033210- 


-275612 


(350^m) 




03 


32 


11, 


.27 


-27 


56 


12.1 


42 


± 


10 


3. 


,96 


273 


BLAST 


J033234- 


-280043 


(350nm) 




03 


32 


34, 


.46 


-28 


00 


36.3 


42 




10 


3. 


.94 


276 


BLAST 


J033313- 


-273709 


(350Hm) 




03 


33 


13, 


.26 


-27 


37 


14.9 


43 


± 


10 


3. 


.93 


291 


BLAST 


J033136- 


-274933 


(350Hm) 




03 


31 


36 


.08 


-27 


49 


40.9 


41 


± 


10 


3. 


,89 


298 


BLAST 


J033312- 


-275611 


(350^m) 


106,107 


03 


33 


13, 


.37 


-27 


56 


04.2 


41 


± 


10 


3, 


,88 


317 












03 


32 


29, 


.08 


-28 


01 


53.8 


41 


± 


10 


3, 


,84 


322 


BLAST 


J033148- 


-280212 


(350Hm) 


29 


03 


31 


48, 


.41 


-28 


02 


20.9 


41 


± 


10 


3, 


,83 


340 


BLAST 


J033132- 


-274314 


(350Hm) 




03 


31 


32, 


.37 


-27 


43 


16.5 


41 


± 


10 


3. 


,80 


362 












03 


33 


04, 


.96 


-27 


37 


26.0 


40 


± 


10 


3. 


,76 
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Table C3. 500 |^m peak list produced from matched-filtered maps. Columns have the same meaning as in Table [CT] On average, each 
of these peaks is a blend of 3.7 sources from our matching catalogue. 



SID Previously Pubhshed Match R.A. Dec. Flux Density SNR 



500 




Name 




ID 




m 




(° 




") 


(mjy) 






5 


BLAST 


J033328- 


275659 


(500n.m) 


118 


03 


33 


28 


17 


—27 


56 


57 





62 


± 


8 


7 


43 


7 


BLAST 


J033311- 


275610 


(500n.m) 


106 


03 


33 


11 


77 


—27 


56 


11 


2 


56 


± 


8 


6 


81 


8 


BLAST 


J033207- 


275811 


(500n.m) 


44 


03 


32 


08 


12 


—27 


58 


05 





53 


± 


8 


6 


51 


9 


BLAST 


J033129- 


275545 


(500n.m) 


6,7,8 


03 


31 


29 


24 


—27 


55 


49 





53 


± 


8 


6 


41 


11 


BLAST 


J033128- 


275713 


(500n.m) 


10,11 


03 


31 


29 


51 


-27 


57 


22 





53 


± 


8 


6 


28 


12 


T AO rn 

BLAST 


J033258- 


274325 


(500n.m) 


92 


03 


32 


58 


68 


—27 


43 


27 


9 


52 


± 


8 


6 


26 


13 


BLAST 


J033215- 


275027 


(500n.m) 




03 


32 


15 


44 


—27 


50 


30 


6 


49 


± 


8 


6 


08 


14 


T AO rn 

BLAST 


J033153- 


273921 


(500n.m) 


34,35 


03 


31 


52 


05 


—27 


39 


21 


8 


50 


± 


8 


6 


08 


17 


BLAST 


J033128- 


273942 


(500n.m) 


2,3,4 


03 


31 


28 


18 


-27 


39 


42 


3 


53 


± 


8 


6 


01 


20 


BLAST 


J033256- 


280102 


(500n.m) 


91 


03 


32 


56 


40 


-28 


01 


05 


4 


48 


± 


8 


5 


84 


21 


BLAST 


J033217- 


275212 


(500n.m) 


55,56 


03 


32 


18 


29 


-27 


52 


12 


7 


47 




8 


5 


74 


23 


BLAST 


J033321- 


275510 


(500n.m) 


117 


03 


33 


22 


24 


-27 


55 


14 


4 


47 




8 


5 


70 


24 


BLAST 


J033318- 


274926 


(500n.m) 


113,114,115 


03 


33 


17 


67 


-27 


49 


26 


6 


46 


± 


8 


5 


67 


26 










70,71 


03 


32 


37 


60 


—27 


35 


31 


5 


46 


± 


8 


5 


61 


31 


BLAST 


J033220- 


275631 


(500n.m) 


57 


03 


32 


20 


86 


—27 


56 


26 


3 


43 


± 


8 


5 


27 


33 


BLAST 


J033212- 


275558 


(500n.m) 


50 


03 


32 


13 


11 


-27 


56 


03 


4 


43 


± 


8 


5 


24 


35 


T AO rn 

BLAST 


J033156- 


274511 


(500n.m) 




03 


31 


58 


33 


—27 


45 


13 


5 


42 


± 


8 


5 


20 


37 


BLAST 


J033249- 


274227 


(500n.m) 


83,84,85 


03 


32 


47 


94 


-27 


42 


25 


3 


43 


± 


8 


5 


14 


41 


BLAST 


J033253- 


273759 


(500n.m) 


90 


03 


32 


52 


93 


-27 


37 


59 


7 


41 




8 


5 


05 


43 


BLAST 


J033301- 


275625 


(500n.m) 




03 


33 


00 


48 


-27 


56 


08 


5 


41 




8 


5 


05 


44 


BLAST 


J033309- 


275125 


(500iJ.m) 


104 


03 


33 


09 


35 


-27 


51 


20 


1 


41 


± 


8 


5 


05 


48 


BLAST 


J033135- 


275448 


(500n.m) 


13,14,15,16 


03 


31 


36 


29 


-27 


54 


42 


9 


41 


± 


8 


4 


99 


53 


T~i T AO rn 

BLAST 


J033229- 


274426 


(500n.m) 


61,62 


03 


32 


29 


05 


—27 


44 


30 


6 


41 


± 


8 


4 


96 


56 


BLAST 


J033235- 


275518 


(500n.m) 


66 


03 


32 


35 


23 


-27 


55 


14 


6 


41 




8 


4 


91 


57 










108 


03 


33 


14 


92 


-27 


45 


32 





40 


± 


8 


4 


90 


58 


BLAST 


J033210- 


273725 


(500n.m) 


49 


03 


32 


10 


34 


-27 


37 


20 





40 


± 


8 


4 


87 


64 


BLAST 


J033129- 


275918 


(500n.m) 


9 


03 


31 


29 


49 


—27 


59 


14 


7 


41 




8 


4 


77 


68 


BLAST 


J033211- 


275210 


(500n.m) 


46 


03 


32 


10 


34 


—27 


51 


54 





39 




8 


4 


74 


84 


BLAST 


J033253- 


274459 


(500n.m) 




03 


32 


54 


41 


-27 


44 


52 


3 


37 


± 


8 


4 


58 


88 


BLAST 


J033153- 


274943 


(500n.m) 




03 


31 


52 


34 


-27 


49 


27 


5 


38 




8 


4 


57 


91 


T AO rn 

BLAST 


J033213- 


274302 


(500n.m) 


48 


03 


32 


13 


32 


—27 


43 


01 


3 


38 




8 


4 


55 


103 












03 


32 


34 


40 


-27 


50 


17 


3 


36 


± 


8 


4 


43 


104 


BLAST 


J033229- 


274314 


(500n.m) 


63 


03 


32 


29 


47 


-27 


43 


08 


2 


36 


± 


8 


4 


42 


109 


BLAST 


J033256- 


274539 


(500n.m) 




03 


32 


56 


26 


-27 


45 


49 


6 


35 


± 


8 


4 


37 


128 


T~i T AO rn 

BLAST 


J033318- 


274115 


(500n.m) 




03 


33 


19 


16 


—27 


41 


15 


3 


35 


± 


8 


4 


28 


133 


BLAST 


J033238- 


275638 


(500n.m) 


72 


03 


32 


38 


87 


-27 


56 


44 


9 


35 


± 


8 


4 


26 


137 


BLAST 


J033323- 


274900 


(500n.m) 




03 


33 


23 


47 


-27 


48 


53 


7 


35 


± 


8 


4 


26 


144 










86 


03 


32 


49 


97 


-27 


58 


33 


2 


35 


± 


8 


4 


24 


156 


BLAST 


J033137- 


273743 


(500n.m) 




03 


31 


37 


89 


—27 


37 


30 


7 


36 




8 


4 


16 


157 


BLAST 


J033157- 


275930 


(500n.m) 


40 


03 


31 


56 


89 


—27 


59 


26 


3 


34 




8 


4 


15 


162 


T AO rn 

BLAST 


J033143- 


274817 


(500n.m) 


22 


03 


31 


43 


92 


—27 


48 


19 


2 


33 


± 


8 


4 


14 


176 


BLAST 


J033300- 


274852 


(500n.m) 


94 


03 


32 


59 


71 


—27 


48 


53 


8 


32 


± 


8 


4 


06 


1 '7'7 










zy 


HQ 

Uo 


31 


48 


o4 


OQ 
— ZO 


no 


03 


3 


34 


± 


8 


A 

4 


UO 


191 












03 


32 


13 


32 


-27 


59 


08 


4 


33 


± 


8 


4 


03 


199 


BLAST 


J033247- 


275415 


(500n.m) 


82 


03 


32 


47 


54 


-27 


54 


14 





33 


± 


8 


4 


01 


200 










98 


03 


33 


06 


56 


-28 


00 


46 


8 


33 


± 


8 


4 


01 


202 


BLAST 


J033154- 


275343 


(500n.m) 


38,39 


03 


31 


54 


94 


-27 


53 


35 


8 


33 




8 


4 


00 


203 


BLAST 


J033215- 


275901 


(500n.m) 


52 


03 


32 


17 


57 


-27 


58 


46 


5 


33 


± 


8 


4 


00 


205 


BLAST 


J033153- 


274943 


(500n.m) 


36 


03 


31 


54 


51 


-27 


49 


43 


2 


33 


± 


8 


3 


99 


210 










42 


03 


32 


06 


41 


-28 


01 


07 


1 


32 


± 


8 


3 


97 


222 


BLAST 


J033309- 


275347 


(500n.m) 


97 


03 


33 


09 


20 


-27 


54 


00 


5 


32 


± 


8 


3 


94 


223 


BLAST 


J033222- 


274555 


(500n.m) 




03 


32 


22 


75 


-27 


45 


52 


9 


31 


± 


8 


3 


93 


248 


BLAST 


J033145- 


274148 


(500n.m) 




03 


31 


46 


11 


-27 


41 


44 


2 


31 


± 


8 


3 


87 


250 


BLAST 


J033235- 


274931 


(500n.m) 


67 


03 


32 


35 


56 


-27 


49 


04 


7 


32 


± 


8 


3 


87 


253 


BLAST 


J033307- 


274001 


(500n.m) 




03 


33 


07 


81 


-27 


39 


36 


2 


32 


± 


8 


3 


86 


261 


BLAST 


J033148- 


280315 


(500n.m) 




03 


31 


49 


06 


-28 


03 


13 


2 


31 




8 


3 


84 


267 


BLAST 


J033243- 


273914 


(500n.m) 


78 


03 


32 


43 


67 


-27 


39 


38 


1 


31 


± 


8 


3 


83 
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Table C4. 870 ^im peaks from [Weiss et aL] ( |2009^ that land within the coverage of the matching catalogue. The columns have the same 
moaning as in Table |C1| although the SNR is measured with respect to the combined instrumental and confusion noise. On average, 
each of these peaks is a blend of 2.2 sources in our matching catalogue. In addition, the 33 peaks marked with an asterisk in the 'Match 
ID' column have no matches in any of the BLAST peak catalogues (Tables [ClJC3]l . These peaks are used to measure priors, but are not 
otherwise analyzed in this paper. 



SID Previously Published Match R.A. Dec. Flux Density SNR 



870 




Name 




ID 




m 


= ) 


{° 




") 


(mJy) 






1 


LESS 


J033314.3- 


■275611 


107 


03 


33 


14, 


.26 


-27 


56 


11.2 


14.7 ± 1.2 


12.48 


2 


LESS 


J033302.5- 


■275643 


95 


03 


33 


02, 


.50 


-27 


56 


43.6 


12.2 ± 1.2 


10.31 


3 


LESS 


J033321.5- 


■275520 


117 


03 


33 


21, 


.51 


-27 


55 


20.2 


11.9 ± 1.2 


10.12 


4 


LESS 


J033136.0- 


■275439 


13,14,15 


03 


31 


36, 


.01 


-27 


54 


39.2 


11.2 ± 1.1 


9.72 


5 


LESS 


J033129.5- 


■275907 


9 


03 


31 


29, 


.46 


-27 


59 


07.3 


10.1 ± 1.2 


8.45 


6 


LESS 


J033257.1- 


■280102 


91 


03 


32 


57, 


.14 


-28 


01 


02.1 


9.8 




X 


2 


8.16 


7 


LESS 


J033315.6- 


■274523 


108 


03 


33 


15, 


.55 


-27 


45 


23.6 


9.4 


± 


1 


,2 


7.89 


9 


LESS 


J033211.3- 


■275210 


46 


03 


32 


11, 


.29 


-27 


52 


10.4 


9.4 




X 


2 


7.71 


10 


LESS 


J033219.0- 


■275219 


56 


03 


32 


19, 


.02 


-27 


52 


19.4 


9.3 


± 


1 


2 


7.60 


11 


LESS 


J033213.6- 


■275602 


50 


03 


32 


13, 


.58 


-27 


56 


02.5 


9.2 




X 


2 


7.58 


12 


LESS 


J033248.1- 


■275414 


82 


03 


32 


48, 


.12 


-27 


54 


14.7 


8.9 


± 


1 


,2 


7.24 


13 


LESS 


J033249.2- 


■274246 


84,85 


03 


32 


49, 


.23 


-27 


42 


46.6 


8.9 




X 


2 


7.21 


14 


LESS 


J033152.6- 


■280320 


33 


03 


31 


52, 


.64 


-28 


03 


20.4 


9.5 


± 


1 


.3 


7.19 


16 


LESS 


J033218.9- 


■273738 




03 


32 


18, 


.89 


-27 


37 


38.7 


8.2 


-|- 


X 


2 


6.87 


17 


LESS 


J033207.6- 


■275123 




03 


32 


07, 


.59 


-27 


51 


23.0 


7.8 


-I- 


1 


,2 


6.36 


18 


LESS 


J033205.1- 


■274652 


43 


03 


32 


05, 


.12 


-27 


46 


52.1 


7.7 


-|- 


X 


2 


6.30 


19 


LESS 


J033208.1- 


■275818 


44 


03 


32 


08, 


.10 


-27 


58 


18.7 


7.5 


± 


1 


2 


6.19 


20 


LESS 


J033316.6- 


■280018 




03 


33 


16, 


.56 


-28 


00 


18.8 


7.5 


-|- 


X 


2 


6.15 


25 


LESS 


J033157.1- 


■275940 


40 


03 


31 


57, 


.05 


-27 


59 


40.8 


7.0 


-I- 


1 


,2 


5.83 


26 


LESS 


J033136.9- 


■275456 


15,16 


03 


31 


36, 


.90 


-27 


54 


56.1 


6.8 


-|- 


X 


2 


5.80 


28 


LESS 


J033302.9- 


■274432 




03 


33 


02, 


.92 


-27 


44 


32.6 


7.0 


-|- 


X 


2 


5.62 


31 


LESS 


J033150.0- 


■275743 




03 


31 


49, 


.96 


-27 


57 


43.9 


6.7 


-|- 


X 


2 


5.51 


32 


LESS 


J033243.6- 


■274644 




03 


32 


43, 


.57 


-27 


46 


44.0 


6.8 


-|- 


X 


2 


5.46 


33 


LESS 


J033149.8- 


■275332 




03 


31 


49, 


.78 


-27 


53 


32.9 


6.8 


-|- 


X 


2 


5.46 


34 


LESS 


J033217.6- 


■275230 


55 


03 


32 


17, 


.64 


-27 


52 


30.3 


6.8 


-|- 


X 


2 


5.40 


36 


LESS 


J033149.2- 


■280208 


29 


03 


31 


49, 


.15 


-28 


02 


08.7 


o.y 


_i_ 
m 


1 


.O 


5.36 


38 


LESS 


J033310.2- 


■275641 




03 


33 


10, 


.20 


-27 


56 


41.5 


6.4 


± 


1, 


.2 


5.22 


40 


LESS 


J033246.7- 


■275120 




03 


32 


46, 


.74 


-27 


51 


20.9 


6.4 


± 


1, 


.2 


5.18 


42 


LESS 


J033231.0- 


■275858 


65 


03 


32 


31, 


.02 


-27 


58 


58.1 


6.4 


± 


1, 


.2 


5.13 


43 


LESS 


J033307.0- 


■274801 


99,100 


03 


33 


07, 


.00 


-27 


48 


01.0 


6.4 


± 


1, 


.3 


5.12 


45 


LESS 


J033225.7- 


■275228 




03 


32 


25, 


.71 


-27 


52 


28.5 


6.3 


± 


1, 


.2 


5.10 


49 


LESS 


J033124.5- 


■275040 




03 


31 


24, 


.45 


-27 


50 


40.9 


6.6 


± 


1, 


.3 


5.05 


50 


LESS 


J033141.2- 


■274441 


19,20 


03 


31 


41, 


.15 


-27 


44 


41.5 


6.1 


± 


1, 


.2 


5.02 


51 


LESS 


J033144.8- 


■274425 




03 


31 


44, 


.81 


-27 


44 


25.1 


6.2 


± 


1, 


.2 


5.01 


52 


LESS 


J033128.5- 


■275601 


7 


03 


31 


28, 


.51 


-27 


56 


01.3 


6.2 


± 


1, 


.2 


4.94 


53 


LESS 


J033159.1- 


■275435 




03 


31 


59, 


.12 


-27 


54 


35.5 


6.2 


± 


1, 


.2 


4.93 


55 


LESS 


J033302.2- 


■274033 




03 


33 


02, 


.20 


-27 


40 


33.6 


6.1 


± 


1, 


.2 


4.90 


56 


LESS 


J033153.2- 


■273936 


35 


03 


31 


53, 


.17 


-27 


39 


36.1 


6.0 


± 


1, 


.2 


4.89 


57 


LESS 


J033152.0- 


■275329 




03 


31 


51, 


.97 


-27 


53 


29.7 


6.1 


± 


1, 


.3 


4.87 


59 


LESS 


J033303.9- 


■274412 




03 


33 


03, 


.87 


-27 


44 


12.2 


6.0 


± 


1, 


.2 


4.77 


60 


LESS 


J033317.5- 


■275121 




03 


33 


17, 


.47 


-27 


51 


21.5 


5.8 


± 


1, 


.2 


4.75 


61 


LESS 


J033245.6- 


■280025 




03 


32 


45, 


.63 


-28 


00 


25.3 


5.9 


± 


1, 


.2 


4.73 


63 


LESS 


J033308.5- 


■280044 


98 


03 


33 


08, 


.46 


-28 


00 


44.3 


6.0 


± 


1, 


.3 


4.71 


64 


LESS 


J033201.0- 


■280025 




03 


32 


01, 


.00 


-28 


00 


25.6 


5.8 


± 


1, 


.2 


4.70 


65 


LESS 


J033252.4- 


■273527 


* 


03 


32 


52, 


.40 


-27 


35 


27.7 


5.9 


± 


1, 


.3 


4.67 


67 


LESS 


J033243.3- 


■275517 


76,77 


03 


32 


43, 


.28 


-27 


55 


17.9 


5.9 


± 


1, 


.3 


4.67 


68 


LESS 


J033233.4- 


■273918 




03 


32 


33, 


.44 


-27 


39 


18.5 


5.8 


± 


1, 


.2 


4.65 


69 


LESS 


J033134.3- 


■275934 




03 


31 


34, 


.26 


-27 


59 


34.3 


5.7 




1, 


.2 


4.65 


70 


LESS 


J033144.0- 


■273832 




03 


31 


43, 


.97 


-27 


38 


32.5 


5.7 


± 


1, 


.2 


4.64 


72 


LESS 


J033240.4- 


■273802 


73 


03 


32 


40, 


.40 


-27 


38 


02.5 


5.7 


± 


1, 


.2 


4.63 
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SID Previously Published Match R.A. Dec. Flux Density SNR 



870 




Name 




ID 




m 


(° 




") 




(mjy) 






73 


LESS 


J033229.3- 


-275619 




03 


32 


29.33 


—27 


56 


19.3 


5. 


.8 


± 


1, 


.2 


4.63 


74 


LESS 


J 033309. 3- 


-274809 


101,102,103 


03 


33 


09.34 


—27 


48 


09.9 


5. 


.8 


± 


1, 


.3 


4.62 


75 


LESS 


J033126.8- 


-275554 


6 


03 


31 


26.83 


—27 


55 


54.6 


5, 


.8 


± 


1 


.3 


4.61 


77 


LESS 


J033157.2- 


275633 




03 


31 


57.23 


—27 


56 


33.2 


5, 


,5 


± 


1 


.2 


4.42 


79 


LESS 


J033221.3- 


-275623 


57 


03 


32 


21.25 


—27 


56 


23.5 


5, 


,5 


± 


1 


.2 


4.40 


80 


LESS 


J033142.2- 


-274834 


21,22 


03 


31 


42.23 


—27 


48 


34.4 


5, 


,4 


± 


1 


.2 


4.38 


81 


LESS 


J033127.5- 


-274440 


1 


03 


31 


27.45 


—27 


44 


40.4 


5, 


,7 


± 


1 


.3 


4.38 


82 


LESS 


J033253.8- 


-273810 


90 


03 


32 


53.77 


—27 


38 


10.9 


5, 


,3 


± 


1 


.2 


4.35 


84 


LESS 


J033154.2- 


-275109 




03 


31 


54.22 


—27 


51 


09.8 


5. 


.5 


± 


1, 


.3 


4.33 


88 


LESS 


J033155.2- 


-275345 


38,39 


03 


31 


55.19 


—27 


53 


45.3 


5. 


.4 




1, 


.3 


4.28 


89 


LESS 


J033248.4- 


-280023 




03 


32 


48.44 


—28 


00 


23.8 


5. 


.3 




1, 


.2 


4.25 


90 


LESS 


J033243.7- 


-273554 




03 


32 


43.65 


—27 


35 


54.1 


5, 


.4 


± 


1 


.3 


4.23 


91 


LESS 


J033135.3- 


274033 




03 


31 


35.25 
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Table C5. Matches between the external catalogue with each of the 4 submm peak lists. The first column gives a unique short ID 
that is used throughout this paper. Since all radio sources and the FIDEL catalogue use the SIMPLE survey {Spitzer IRAC/MUSYC 
Public Legacy in ECDF-S; [Gawiser et al.'2006') as a positional prior, the SIMPLE IDs for each source are given in column two. The only 
exceptions are 870 |^m peaks with no identification in the matching catalogue; these peaks have simply been added to the catalogue with 
the SIMPLE ID replaced by the 870 |J.m SID referring to Table |C4| (in boldface brackets). The third and fourth columns give the IRAC 
J2000 coordinates of each source (or the 870 |xm position when unavailable). There are then 4 blocks of 5 columns for each submm band 
giving: the submm SID for each respective submm peak list to which the source was matched (Tables |C1]C4| |; the radial separation in 
arcsec; the likelihood ratio; and finally the 24 |j.m and 1.4 GHz P values when available. 
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17 


8 
















90 


6.5 


2.9 





102 





106 






























90 


(82) 


03 


32 


53 


77 


-27 


38 


10 


9 
















52 


31.8 












41 


15.8 












5 


0.0 
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ID SIMPLE R.A. Dec. 250 ^m 350 ^m 500 nm 870 ^m 

SID Ar LR P24 Pr SID Ar LR P24 P, SID Ar LR P24 Pr SID Ar LR P24 Pr 
(h m S) (O , //) (//) (//) (//) 



91 


(8) 


03 


32 


57.14 


— 28 


01 


02 


1 


































20 


10.3 






























92 


38977 


03 


32 


57.63 


— 27 


43 


18 





27 


11.4 


2. 


4 





100 





136 


















































93 


50259 


03 


32 


59.32 


— 27 


35 


34 


1 


130 


5.3 


21 








030 





019 


















































94 


29432 


03 


32 


59.32 


— 27 


48 


58 


6 


































176 


7.0 


10.2 





029 





094 




















95 


17004 


03 


33 


01.59 


— 27 


56 


49 


4 


















64 


9.9 


2 


5 





240 





142 


































96 


37597 


03 


33 


06.16 


— 27 


44 


15 


3 


















115 


1.8 


1 


3 





Oil 





016 


































97 


21942 


03 


33 


07.78 


-27 


53 


51 


3 


































222 


20.9 


1.4 





226 





335 




















98 
99 


(63) 
31482 


03 
03 


33 
33 


08.46 
06.63 


-28 
-27 


00 
48 


44 
02 


3 



































200 


25.2 












43 


5 





2 


6 





099 





396 


100 


30926 


03 


33 


07.25 


-27 


48 


08 


4 
















































43 


8 


1 


1 


1 





074 





189 


101 


31537 


03 


33 


09.15 


-27 


48 


16 


8 


















145 


8.7 


1 


4 





278 





169 
















74 


7 


4 


4 


9 





222 





126 


102 


30714 


03 


33 


09.71 


-27 


48 


01 


6 


52 


14.7 


27 


5 





031 





052 


145 


22.0 


2 


6 
























74 


9 


7 


11 


.3 





015 





025 


103 


30905 


03 


33 


09.77 


-27 


48 


21 



















































74 


12 


.5 


2 


2 










104 


26260 


03 


33 


10.12 


-27 


51 


24 


8 


































44 


11.3 


1.3 





263 





237 




















105 


41860 


03 


33 


11.79 


-27 


41 


38 


4 


















31 


10.2 


6 


5 





256 





035 


































106 


18562 


03 


33 


12.62 


-27 


55 


51 


7 


















298 


15.9 


2 


1 
























96 


7 


5 


1 


2 





024 





073 


107 


(1) 


03 


33 


14.26 


-27 


56 


11 


2 


















298 


13.8 














7 


32.9 












1 




















108 


35704 


03 


33 


15.42 


-27 


45 


24 





















98 


8.8 


2 


4 





113 





113 


57 


10.3 


1.6 





142 





144 


7 


1 


8 


4 


4 





009 





009 


109 


26283 


03 


33 


16.51 


-27 


50 


39 


4 


72 


4.2 


39 


5 





002 





007 
































108 


6 


4 


4 


2 





004 





015 


110 


41704 


03 


33 


16.93 


-27 


41 


21 


4 


185 


10.7 


27 


3 





017 





039 


















































111 


34567 


03 


33 


17.78 


-27 


46 


05 


9 


54 


2.6 


7. 







012 





Oil 


43 


11.1 


5 


9 





122 





115 


































112 


34104 


03 


33 


17.78 


-27 


46 


23 


6 


















43 


12.9 


1 


8 





294 





150 


































113 


28590 


03 


33 


17.43 


-27 


49 


48 


7 


































24 


22.3 


2.7 





300 





177 




















114 


29761 


03 


33 


17.81 


-27 


49 


10 


4 


































24 


16.3 


1.1 





538 





375 




















115 


28465 


03 


33 


18.69 


-27 


49 


40 


1 


































24 


19.2 


2.0 





549 





128 




















116 


20681 


03 


33 


18.90 


-27 


54 


33 


5 


67 


9.9 


1. 







129 





141 


















































117 


19773 


03 


33 


21.48 


-27 


55 


20 


3 
















































3 





4 


2 


4 





010 





015 


118 


17028 


03 


33 


28.56 


-27 


56 


54 


1 


















14 


15.2 


1 








273 





313 


5 


5.9 


5.0 





075 





093 


120 


2 


3 


11 


.2 





015 





019 
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Table C6. Maximum-likelihood submm photometry, SED model fits and redshifts. Uncertainties in flux densities have the estimated 
confusion noise (Tables [T]l added in quadrature to the instrumental noise. The temperatures, Tq^s and 10-1000 |J.m total infrared (TIR) 
fluxes. Stir are for observed-frame modified blackbody fits of the form Sv oc i'^'^ B,j(Ta\^s) (rest-frame values shown in Fig. \El\ . 
The fractional uncertainties in the TIR fluxes, AStir/Stir, are produced from the same Monte Carlo simulation used to measure 
uncertainties in the temperatures. SED fits are not provided for confused sources (see SEDs with 'C indicated in Fig. When 
redshifts are available, the TIR luminosity, LtiRi is also calculated, but integrating in the rest-frame, and using a more realistic SED fit 
from the library of [Dale et al.| |2001| in order to estimate emission in the mid/far-IR (fits shown in Fig. |E1| |. Redshifts in boldface are 
optical spectroscopic measurements, redshifts in regular face are optical photometric estimates, and redshifts in brackets are IRAC-based 
photometric estimates. The superscripts indicate where the redshifts were found: 'P' from the composite catalogue of |Pascale et alT] 
( |2009| l; T from |Ivison et al.| l |2Mol l; 'D' from |Dunlop et al.|poTol l; and 'C from |Casey et al.|p010l l. 



ID 




■Saso 


Ssoo 


5870 




Stir 


AStir/Stir 


^TIR 


Rcdshift 




(mjy) 


(mJy) 


(mJy) 


(mJy) 


(K)' 


(logio Wm-2) 




(log^Q Lq) 


{^) 


1 


38±18 


41±15 


9±13 


5.1±1.2 


lO.Sitl.l 


-15.3 


0.36 


12.4 


(1.9)P 


2 


17±33 


-40±30 


-92±31 


1.2±2.2 










l.lP 


3 


40±31 


95±29 


131±30 


2.3±2.1 










0.8^ 


4 


26±32 


-2±30 


0±31 


2.1±1.5 










I.7P 


5 


62±30 


56±29 


29±30 


2.2±1.4 










2.6^ 


6 


13±21 


30±19 


42±18 


3.6±1.2 


9.7±1.0 


-15.5 


0.35 


12.6 


2.7P 


7 

8 


-20±34 
93±30 


-7±33 
29±29 


-28±39 
49±31 


4.1±1.3 
-0.2±1.3 










(1.9)P 
O.677I 


9 


69±18 


47±15 


51±13 


9.6±1.2 


9.9±0.7 


-15.1 


0.24 






10 


109±24 


85±20 


60±18 


3.3±1.2 










1.482^ 


11 


-14±24 


-25±20 


-14±18 


2.8±1.2 










O.4P 


12 


68±18 


40±15 


22±13 


-0.4±1.2 


17.7±0.5 


-14.8 


0.15 


12.4 


I.2P 


13 


-39±47 


-19±52 


29±50 


4.1±1.8 












14 


31±85 


57±101 


-5±101 


2.9±2.5 










(1.2)P 


15 


16±58 


20±68 


18±67 


6.3±1.8 










O.lP 


16 


24±22 


-8±20 


14±19 


4.4±1.2 


6.1±1.9 


-16.2 


0.72 


11.6 


2.1P 


17 


41±18 


43±15 


24±13 


4.2±1.2 


11.4±1.8 


-15.2 


0.48 


12.6 


2.O24P 


18 


55±18 


20±15 


13±13 


3.6±1.2 


13.0±4.5 


-15.1 


0.90 


13.1 


3.0P 


19 


32±29 


51±26 


51±23 


2.5±1.3 










(0.6)P 


20 


28±29 


-2±26 


-23±23 


3.9±1.3 










O.4P 


21 


21±27 


13±23 


24±22 


2.9±1.3 










I.2P 


22 


-9±27 


4±23 


3±22 


3.2±1.3 










(2.2)P 


23 


21±19 


37±16 


28±15 


2.5±1.2 


11.4±2.1 


-15.4 


0.67 


12.1 


(1.5)P 


24 


28±19 


25±16 


0±15 


1.7±1.2 


21.6±1.3 


-14.7 


0.19 


9.9 


O.O95P 


25 


65±18 


33±15 


30±13 


1.4±1.2 


15.7±1.0 


-14.9 


0.24 


12.2 


Ll^ 


26 


77±18 


22±15 


22±13 


3.9±1.2 


14.4±3.7 


-14.9 


0.90 






27 


78±18 


45±15 


22±13 


0.4±1.2 


23.2±0.9 


-14.3 


0.09 


11.6 


0.364^ 


28 


40±18 


38±15 


13±13 


3.6±1.2 


11.9±1.8 


-15.2 


0.63 


11.3 


0.6^ 


29 


35±18 


39±15 


30±13 


6.5±1.2 


9.6±1.1 


-15.3 


0.41 


12.5 


(2.1)^ 


30 


84±18 


53±15 


33±13 


3.2±1.2 


18.2±0.4 


-14.5 


0.08 


12.0 


0.620^ 


31 


18±22 


-11±20 


25±19 


1.7±1.2 


9.2±2.9 


-16.0 


0.90 


12.2 


(3.0)P 


32 


57±23 


79±20 


15±19 


3.4±1.2 


13.5±1.8 


-14.9 


0.42 


12.6 


I.6O5C 


33 


29±18 


33±15 


22±13 


8.1±1.3 


8.4±1.0 


-15.4 


0.33 


11.7 


(1.2)P 


34 


42±19 


29±17 


25±15 


1.1±1.2 


14.5±1.4 


-15.1 


0.45 


12.8 


2.342'^ 


35 


74±20 


50±17 


38±15 


5.2±1.2 


12.3±1.2 


-14.9 


0.37 


13.2 


2.9P 


36 


60±18 


60±15 


47±13 


3.2±1.2 










(2.1)P 


37 


54±18 


43±15 


24±13 


2.6±1.2 










0.6P 


38 


30±27 


35±24 


49±22 


1.7±1.3 












39 


14±27 


5±24 


-15±22 


3.4±1.3 










1.8P 


40 


11±18 


31±15 


26±13 


6.1±1.2 


8.3±0.9 


-15.6 


0.33 


11.9 


(1.6)P 


41 


17±18 


31±15 


0±13 


1.3±1.2 


18.5±1.4 


-15.0 


0.30 


9.9 


O.lP 


42 


30±18 


12±15 


11±13 


2. Oil. 2 


12.8±6.0 


-15.4 


0.90 


10.1 


O.2P 


43 


68±18 


58±15 


20±13 


6.3±1.2 


11.7±1.0 


-15.0 


0.32 


12.9 


2.252^ 


44 


55±18 


60±15 


65±13 


6.0±1.2 


10.8±0.7 


-15.1 


0.29 


9.7 


(O.l)P 


45 


37±18 


5±15 


5±13 


-0.1±1.2 


32.6±7.8 


-14.1 


0.59 


12.5 


O.7P 


46 


22±18 


55±15 


43±13 


8.5±1.2 


8.8±0.6 


-15.3 


0.32 


12.2 


(1.6)P 


47 


42±18 


-15±16 


-10±14 


1.7±1.2 


17.5±1.3 


-15.2 


0.24 


11.2 


0.6P 


48 


15±18 


38±16 


43±14 


2. Oil. 2 


ll.litS.O 


-15.4 


0.90 


12.1 


1.7° 


49 


30±18 


52±15 


32±13 


2.1±1.2 


12.5±1.4 


-15.2 


0.53 


12.3 


1.565P 


50 


16±18 


10±15 


35±13 


8.1±1.2 


7.1±1.3 


-15.7 


0.52 


10.7 


(0.6)P 
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ID 


5*250 


5350 


^500 


'S'sTO 


'^obs 


5'tir 


^5tir,/5tir 


^TIR 


Redshift 




(mJy) 


(mJy) 


(mJy) 


(mJy) 


(K) 


(logi„ Wm-2) 




(logio L©) 


(2) 


51 


48dil8 


55dil5 


12dil3 


-0.6dil.2 


17.0i0.5 


-14.9 


0.21 


12.4 


1.324^ 


52 


33dil9 


33dil6 


21dil4 


0.2dil.2 


18.2i5.1 


-14.9 


0.90 


11.6 


0.6^ 


53 


74dil9 


51dil6 


17dil4 


2.1dil.2 


21.0i0.7 


-14.4 


0.11 


10.5 


0.1? 


54 


56dil8 


26dil5 


17dil3 


1.5dil.2 


90 <^-\-C\ Q 
zu.o^u . y 


-14.6 


0.15 


10.3 


0.124^ 


55 


46dil9 


24dil6 


30dil4 


3. Oil. 2 


12.5i2.1 


-15.2 


0.56 


11.9 


1.097^ 


56 


35dil9 


50dil7 


25dil5 


8.2dil 2 


9.2il.l 


-15.3 


0.40 


12.6 


2. 3D 


57 


69dil8 


32dil5 


37dil3 


4.4dil 2 


1 9 fi-UO Q 

-L^ .U— 1— Z.i7 


-15.0 


0.89 


12.9 


2.277^ 


oo 


f^i -1-1 >^ 




-O^ lO 


n •^-1-1 9 

-U.O^ 1 .z 


25 2il 4 


-14 4 


n 14 


in 1 


n iP 

u. 1 


59 


48dil8 


26dil5 


19dil3 


l.Odil.2 


24.4il.2 


-14.4 


0.13 


10.3 


0.1? 


fin 


uy^ -Lo 




Q/l -Ul Q 


1 Q-l-1 9 
1 .O^ 1 .z 


1 PC 7-1-9 1 

lO. ( ^Z. 1 


-14 8 


U.OO 


12 3 


1.1^ 


61 


97di20 


42dil8 


33dil6 


1.9dil.2 


Z'i.O^U . o 


-14.1 


0.08 


10.4 


0.077^ 






-LO 


1 -1-1 fi 

-L^ lU 


n 1 -1-1 9 

U. 1^ 1 .z 


24 8il 1 


14 


n 19 


m 1 

J-U. J. 


0.076^ 


63 


64dil8 


73dil5 


43dil3 


4. Oil. 2 


12.7i2.0 


-14.9 


0.56 


12.4 




fid 




ZU^ -LU 


U^ 1^ 


2 5il 2 










o 1 9^;? 


65 


27dil8 


4dil6 


-10dil4 


4.2il 2 










^^z.o; 


fifi 


J-OO^ lo 


/ O^ -LO 




n 1 -1-1 9 

U. 1^ 1 .z 


22 liO 4 


-14 1 


KJ.UO 


Q 7 
y. / 


n n^sD 
u.woo 


67 


54dil8 


35dil5 


40dil3 


2.7il.2 


13.2il.4 


-15.1 


0.25 


13.1 




fis 
uo 




97_|_i Q 
Z / ^ -LO 


lU 


1 Q-l-1 9 
-1 .O^ 1 .z 










u.u 


69 


26±20 


10±18 


-3±16 


1.6il.2 










(2.2)P 


/ u 


j_uz^zy 


U ( ^zu 


yo^zo 


1 zl-l-1 S 

-1.4:^ 1 .0 












71 


-31di29 


-lldi26 


-53di23 


2.7il.8 










0.2? 


79 


/17_|_1 Q 
^ / ^ -Lo 


9Q_|_l r. 
zy^ -LO 


zo^ lo 


-1 1 .z 


99 Q-l-n s 

ZZ.O^U.O 


-14 6 


n 1 n 


111 

J-X.X 


Q 297? 


73 


30dil8 


35dil5 


21dil3 


4.4il 2 


1 n Q-l-9 n 


-15.3 


0.66 


11.4 


0.830^ 


7A 


fi-l-9Q 
u^zy 


oo^zu 


4:y ^zo 


0. 1^ 1 .o 












75 


43di29 


-42di26 


-42di25 


-l.lil.8 










0.1? 




7Q_(_97 

1 y^z I 


41di24 


24di22 


1 Q-l-1 
1 .O^ 1 .0 










2.123*^ 


77 


-24di27 


13di24 


ldi22 


3.8il 5 










2.8? 


7J^ 


UO^ J-O 


-LO 


91 -1-1 
Zl^ lO 


1 Q-l-1 9 
1 .O^ 1 .z 


1 Q 9-l-f^ 7 


-14 6 


u.yu 


19 1 
±z.x 


O.733P 


79 
ou 


55di27 
1 n_|_97 


33di23 

n-i-9'^ 
u^zo 


g-j-22 
-ldi22 


0.4il.9 

n 4-1-1 Q 
1 .y 










Q 279? 
n 97Q? 


81 


94dil8 


52dil5 


21dil3 


2.3il.2 


21.2i0.7 


-14.3 


0.11 


10.4 


0.104? 






Z4:^ ±0 


Z ( ^ lO 


7 7-1-1 9 
1 . i ^ 1 .z 


Q 0-1-9 1 
y .u^z. 1 


1 4 


u.uy 


12 8 


9 aDf 
z.o 


83 


67di21 


70dil8 


32dil6 


2.4il.2 










0.981? 




-28di21 


QQ_|_1 Q 
-OO^ lO 


-22dil6 


7 Q-l-1 9 
/ .O^ 1 .z 










n q^? 


85 


10di20 


7dil8 


34dil6 


2.8il.2 










0.981? 


ou 


QQ_|_1 Q 

yo^ -Lo 


UO^ -LO 


91 -1-1 
Zl^ lO 


2 3il 2 


1 ft 4.-1-0 S 


-14 6 


n 99 
u.zz 




2.326^ 


87 


10dil8 


-22dil5 


3dil3 


0.7il.2 


29.7ill.3 


-15.6 


0.34 


10.3 


0.4? 


oo 


QO-LI Q 
OO^ -Lo 


-LU^ -LO 


lU^ lo 


1 fi-l-1 9 
1 1 .z 


1 Q ft-l-1 Q 
lO.U^ 1 .o 


-14 9 


0.24 


1 1 fi 


o fi9n? 


89 


-5dil8 


43dil5 


8dil3 


1.8il.2 


10. Oil. 8 


-15.7 


0.41 


11.8 




90 


52±18 


51±15 


42±13 


4.6il.2 


11.6il.8 


-15.1 


0.35 






91 


18±18 


40±15 


39±13 


8.9il.2 


8.2i0.8 


-15.4 


0.30 






92 


90±18 


54±15 


38±13 


2.8il.2 


15.7il.7 


-14.7 


0.33 


13.0 


I.9P 


93 


58±18 


25±15 


7±13 


0.9il.2 


18.4i0.8 


-14.8 


0.17 


10.6 


O.2P 


94 


43±18 


13±15 


29±13 


1.3il.2 


16.0i6.0 


-15.0 


0.90 


11.2 


O.5P 


95 


25±18 


39±15 


32±13 


4.7il.2 


9.9il.4 


-15.4 


0.46 


11.9 


I.3P 


96 


51±18 


52±15 


15±13 


0.2il.2 


19.7i4.8 


-14.6 


0.90 


11.4 


O.4P 


97 


36±18 


9±15 


29±13 


1.9il.2 


12.9i2.2 


-15.3 


0.50 


12.1 


I.5P 


98 


-35±18 


-7±15 


13±13 


5.2il.2 


4.4il.3 


-16.5 


0.66 






99 


-30±29 


48±28 


4±33 


3. Oil. 3 










(3.0)P 


100 


70±31 


-27±30 


35±41 


3. Oil. 3 










O.531P 



'Source 82 has an optical photometric redshift limit 2 > 2.8 from Dunlop et al. (20101. Since this limit is significantly larger 
than its IRAC-based photometric redshift Zi = 1.6 we use the limit as the best estimate for this source. 
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ID 
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APPENDIX E: SPECTRAL ENERGY 
DISTRIBUTIONS 

This Appendix contains observed-frame radio-submm- 
(mid/near)-IR SEDs for all 118 proposed identifications to 



the submm peaks (Fig. El I 
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Figure Dl. 1.8' X 1.8' Postage stamps for each of the fitted sources. Columns from left to right: 250 |xm, 350 |xm, 500 |xm, 870 |xm(SNR 
maps, scaled between -3cr (white) and +5(7 (black), 3.6 |^m, and 1.4 GHz (-2 /xjy (white), +5 /.tjy (black)). Small heavy circles indicate 
the source under consideration. Light smaller circles indicate other nearby sources that have been simultaneously fit in the submm bands. 
Large dashed circles on the submm images indicate input positions obtained from local maxima in the match-filtered images. Sources were 
matched to the external catalogue, and the reported redshifts are those that were used in |Pascale et al.| ( |2009 ) when available (subscript 
's' indicates spectroscopic redshifts, 'p' optical photometric redshifts, and 'i' IRAC-based photometric redshifts). The likelihood ratios 
(LR) of the matches to each submm band are indicated in the top-right corners when available. 
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Figure El. Observed-frame spectral energy distributions. The flux densities at all 118 positions are fit simultaneously in each submm 
map. In this way flux densities in blends are divided up. All matched sources have IRAC photometry from SIMPLE, and most sources 
also have 24 |im (and occasional 70 |xm) flux densities from FIDEL. Most sources also have 1.4 GHz radio measurements. The only 
exceptions are sources for which the BLAST data have only been matched to an 870 ^m peak from LESS, whose position has been used 
to re- measure flux densities in the BLAST bands (no SED information given). Rcdshift labels have the same meaning as in Fig. |El| 
Dotted lines are SED fits to the submm and FIR photometry using the libr ary of [Dale et al. (2001). The 'C indicates 42 sources that are 
confused to the point that the submm photometry is unusable, and only the |Dale et al.| l |2001 ) SEDs are fit (to the radio and mid/near-IR 
photometry). For the remaining 83 sources a modified blackbody with emissivity /3 = 2.0 has been fit, with the maximum-likelihood 
rp.fit- frame, temneratures indicated for the 78 sources with redshifts. 
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